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I^iTHODUCCION 


During the last decade a numDer of Soviet scientists 
have investigated so-called non-linear mechanics, ar.i 
among the most active are certainly to he fo-und Krylcfi 
and Bogoliuboff . An extensive bibliography cf their 
contributions to the subject will be fo’unc at the end. 

A cursory reference to it will quickly discicae the 
fact that in one way or another their work is but 
poorly accessible to the American scientific and tech- 
nical public. The present monograph is essentially a 
very condensed English version of their most extensive 
paper (No. 32 of the Bibliography, in Russian) except 
for the last chapter which is practically a small ex- 
tract of their most mathematical production on the sub- 
ject (No. 1 6 of the Bibliography, in Russian). 

Kryloff and Bogoliuboff consider primarily equa- 
tions of the form ^ 

Sh , S t) 

dt^ . dt 

where t is a small positive quantity and f is a power 

series in t, whose coefficients are polynomials in '' 
dx 

X, sin't, cos t. As a matter of fact, generally f 
contains neither t- nor t. Similar equations are well 
known in astronomy and have been the object of system- 
atic investigation by Llnstedt, Gulden, Llapounoff 
and, above all by Poinpare'.. In a general sense, one 
may say that the same methods are applied by Kryloff and 




X. 


SOME NON-LUTOAR QSCIIXATQRY SY^XMS- 


1 - In ttie present: section we will discuss a Tew 
non-linear oscillatory systems and derive t.'he corres- 
ponding dirferential equations. Tliese equations will 
serve later as illustrations for tb.e metliods c 5 approxi 
mat Ion Introduced in tiie sequel. 

2. We "begin with, some conservative (non-dlssipa- 
tlve ) systems - 

(2.1) Oscillating shaTt . Consider a sinait com- 
posed (ideally) of two revolving masses Joined hy a 
non-linear elastic connection. Let a be the 

moments of* inertia of the revolving masses, and 
their angles of rotation. Let iurther M = 0(0.5-62) he 
the angular momentum of the elastic connection repre- 
sented as a Tunctlon of the angle of rotation 0 — 

The equations of* motion for each of the two masses are 


’’ dt^ 


c(e.,-*e2) = o, 

c(0^-02) == o* 


Hence the equation governing 
(2.2) 


d^0 
_ ^ 


dt" 




the oscillations 

c(e) — o. 



Pig. 1 Pig. 2 


is 
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In this relation the fimctlon c(e) Is usually 
given graphically and may have the most diverse form, 
(2.3) Electrical circuit without resistance . 
Consider an electrical oscillating circuit (Fig. i ) con- 
taining an iron core. Let 5 denote the magnetic fluoc, 
i the line current, C the capacity. We then have 

(2.4) j^ldt = o. 

The relation "between J and i is shown in Pig. 2. With 
sufficient accuracy and within certain limits, one may 
represent this relation analytically, for instance as; 

(a. 5) 1 = Af + 

We have then for $ the differential equation 

( 2 . 6 ) ^ 0 . 

dt^ C 

5. In the exanples of oscillating sya-fllms that we 
have examined so far, we have not taken into considera- 
tion friction which causes dissipation of the oscil- 
lations of the system. 

Generally speaking the laws of mecnanical friction 
have been hut little investigated. In practice, one 
chiefly assumes one of the following three: 

a) The force of friction is proportional to the 
velocity (oscillations in the atmosphere). 

h) The force of friction is proportional to the 
square of the velocity (for oscillations in a liquid). 

c) Coulomb’s law: The force of friction is con- 
stant in magnitude hut depends upon the velocity and 
its direction is opposite the velocity (for example in 
slipping of surfaces upon one another). 

^ ^ ^ ^ free ly oaclllating in the.- atmos - 

^ere . it we assume that friction is proportional to 
the velocity, the equation of oscillations will ^ he 
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(5.2) 


^ ^ ^ sin 6 = C 

dt^ dt 1 


^T-^ere A is a proportionality coefficient^ called fric- 
tion coefficient. 

(5*3) Electrical circuit wltii resistance . ¥e 
suppose that the circuit contains an iron core, an 
ohmic resistance and a capacity (Pig. 5). Let 5 he 
the flux, i the current, R the ohjnic resistance, C the 
canacitv. We will have this time 

^ WAWVv — ^ — I 

1 1 u 

Fig. 3 

g + Ri + ^ i icit =.0 

and hence assuming that -(2.5) holds: 

(3.i^) + R(A+5B|)^) ^ = 0. 

dt^ dt C 


i. Up to the present, we have considered oscil- 
lating systems with or without dissipation (friction). 
Since in practice dissipation is always present in 
3 £>rae form in oscillating systems^ the oscillations 
will fail to die down only if the system contains some 
source of energy which may compensate for the loss of 
energy due to dissipation* This condition may he ful- 
filled in two ways. First, the force acting upon the 
oscillating body (due to its connection with the 
source of energy), may possess a definite periodicity. 
The simplest example of oscillations of this type, said 
to be forced, is found In the vibrations of linear 
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systems subjected to a barmionic disturbance: 

( . 1 ) in^„+ A — + kx = P alnat 

dt^ dt 

where m is the mass, x the displacement, A the dissi- 
pation coefficient, k the spring constant, P the ampli- 
tude of the exterior force, at the frequency of the 
disturbance . 

Second, the source of energy Itself may have no 
specific periodicity but its action upon the oscilla- 
ting body appears to introduce into the system a nega- 
tive dissipation which may compensate for the normal 
positive dissipation caused by the dissipative forces. 
Oscillations of this last type, called auto-oscilla- 
tions, are quite wide-spread and have great importance 
in Physics and Technology. 

To obtain some idea of the manner in which auto- 
oscillations arise, we will examine a system with one 
degree of freedom. 

If the oscillations are of rather small amplitude 
we may write down the customaxy linear equation: 

2 

(4.2) m + A^ + kx =0. 

dt 

As is well known the general solution will be 
-dt 

X = ae co3(ttjt+^) 


where a, p are arbitrary constants. 


6 


2m 





2 


Hence if A>0, then the aniplitude of the small oscil- 
lations will die down according to an exponential 

law. If on the contrary A<0 then the small oscillations 
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will expand and t>ie amplitudes will increase expon- 
entially. 

Since for physical reasons the amplitudes carnot 
increase indefinitely, we must suppose that from a 
certain moment the dissipation coefficient changes its 
sign and becomes positive. This fact may te reflected 
in the differential equation of the oscillations, for 
instance, by replacing the constant coefficient A by 
a variable one: 


dx 2 

A =-A B (|f) 

where A>o, B)>0. We thus obtain a differential equation 
due to Rayleigh: 

2 2 

(i{.. 3 ) m + (-A + B(^) ) ^ + kx = C, 

dt 

This equation shows in particular that the dissi- 
pation is negative for small absolute values of ^ and 
positive when its absolute values are large. 

Thus, small oscillations will expand and large 
oscillations will die down. 

The importance of for self-oscillatory 

systems was already brou^t out by Raylei^ in his 
paper: On maintained vibrations (Phil. Mag. S. 5, 

vol. 15, 1883). 

Another Important equation as regards self-oacil- 
latoiy systems, repeatedly investigated by van der Pol 
and going by his name, is 

(1».4) ^ M + X - 0. 

dt^ , dt 

It may be deduced fj^m (4*5) by making the change of 
variables : 
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and setting - r-r r — t • 

\llan 

5. We will now consider some self -oscillatory 
systems . 

(5.1) Blectronic generator . We refer to Fig. h 
for the various designations of currents (written 1 
with subscripts), voltages (written V, E with sub- 
scripts), etc. 




As we know, however, from the theoi*y of electronic 
Isji^ja, the anode current la a definite function of the 
so called directing potential n = Vg + 



( 5 .^) 1 . 


fi^u) - r(Vg 4- 


where D ia a constant factor, 
lainp . 


one ccn 


ictance of one 



Fig. 5 

In practice D is small relatively to unity. A typical 
curve representing the relation (t.l), the so-called 
characteristic of the lamp is shewn in Pig. 5. 

Suhstituting from (5-1) into (5*5) ana in view of 
(5.2) we find 


(5.5) LC 


d^i. 


dt' 


R dt 


di^ 

f (DE^+(M- LD) — 
dt 


Consider now the following quantiries: 

A = ^ • 


Since the directing potential is 

H DE^ + (M-ED)-^ = + V, 

V will clearly be the variable part of this potential 
induced by the vibrations of the current in the oscil- 
lating circuit, and the constant part induced by a 
source of constant current (for instance a battery). 

In view of this, let us apply to both sides of (5-5) 
the operation 

We then obtain for the unknown V a relationship of 
the form 

(5.6) LC^ + V 4- - (M-IJD)f *(E^ + V)i^ = O, 


dt" 


dt 
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If we choose such that it is .he abscissa of the in- 
flexion of the characteristic in Fig. 5, and neglect 
terms of order >3 in the MacLaurin expansion of 
f(E^ + V), then with siiitable choice of a time unit 
(dimensionless time), ( 5 .C) may be reduced to van der 
Pol ^3 equation ( 4.4), Similar considerations lead to a 
Rayleigh differential equation for 1^. 

(5-7) It is .hardly necessary to observe that if 
a harmonic disturbance is superimposed upon any one of 
the preceding systems ^ there is obtained an equation 
for forced oscillations. Thus we may have in relation 
to a van der Pol system an equation 

( 5 . 8 ) -t ( 1 -x®)— + X = F 3 lno 6 t, 

dt 

and likewise for the other systems. 


II. ELEMFNTAHY THEORY OF THE FIRST APPROXIMATION 

6 . All the -examples discussed In the preceding 
chapter lead to equations of the form 

( 6 . 1 ) ^ + F(x,^,t: 0 . 

dt*^ dt 

We propose to investigate more particularly the so 
called quasi-harmonic case, where there are oscillations 
near the sinusoidal: 

X = a sin (9t+4), 
that is to say when we may write 



TT. T^T.-FMEMARY THE 0 R:S< 


RST APPROXIMATION 


( 6 . 2 ) 


P(x,^,t) = 9^x+£f (x,^,t). 


dx 


where t is £. parameter characterizing the smallness of 
the deviation of F from v^^x. Until further notice we 
assume F, and hence also f, free from the explicit 
variable t. The basic differential equation will thus 
be 

(6.3) ^ V^x .tf(x,a2) = 0,- 

dt^ dt 


and this is the equation which we shall Investigate. 

If we endeavor to solve this equation by the usual 
methods of approximation, notably by the method of 
Poisson, we encounter a. classical difficulty which 
baffled the astronomers of the eighteenth century, 
namely the presence of so-called secular terms, or 
terms of the form txa trigonometric function. In the 
same spirit as the astronomers did in their day, we 
shc.ll endeavor to find methods of approximation which 
yield results free from secular terms. 

In the present chapter we shall describe a very 
intuitive method enabling us *to construct an approxi- 
mate solution which will be free from secular terms. 

7 . We first observe that for e. = o, (6.5) has the 
solution 

(7.1) X = a sin (v^t+^), 

(7.2) ^ = aV cos ("^t + ^), 


where the amplitude a and the phase (|> are constant . For 
convenience the term **frequency" will designate V 
rather than the customary 2 tr>^. 

Consider a, ^ , as new unknown functions of the 
time which are to be determined so that (7.1 ) becomes a 
solution of (6.5). We must have first 

(7-5) ^ ^ sin (vt+<|») + a^ cos (vt++) + aV cos(vt + 4»). 
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1 0 

Hence if we wish to preserve ( 7 - 2 ), or we may say If we 
impose ( 7 . 2 ), then 

(7-^0 ^ sin (Vt + 4‘).+ cos (vt+4)) = o. 

Prom these relations we deduce: 

(7.5) coa (Vt+it)) - Va^ ain (Vt+<t)) 

cit"“ dt dt 

- v^^a sin (Vt + <j>), 

and so finally from ( 6 . 5 ): 

( 7 - 6 ) cos (Vt + <t>) - sin {\>t+<}>) 

= - t f (a sin (V t+<j>) ,aj^cos (vt + <j>)). 

By combining with ( 7 - ^ ) there comes 

( 7 . 7 ) ^ = sin (Vt+4.), avcoa (Vt+<1))7 cos (Vt+<>), 

(7.8) ^ sin (Vt+t), aV cos (Vt+4>jJ aln' (Vt+it)) . 

8. Thus instead of the single differential equa- 
tion of the second order ( 6 . 3 ) in the unknown x, we 
have two differential equations of the first order in 
the two unknowns a, <t> . Notice now that the ri^t hand 
aides of (7-7), (7-8) admit with respect to t the 
period T = Moreover* ^ are proportional to 

the small perameter t, sjd that a, <1> will be slowly ^ 
varying functions of thT tSe^'during the period T, and 
as a first approximation we may, therefore, consider 
them as constant. On the strength of this observation, 
we will indicate at once a simple intuitive method for 
constructing an approximate solution of (7-7>, (7.8). 
For this purpose consider the expressions 
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f(a 3in 6,aV cos cb ) cos 6,f(a sin 6,aP cos <b) 3in 6, 
md let U3 expand them in Fourier series. We find 
'f(a sin 6, a>^ cos cb ) cos 6 ( a )+ a ) 


( 8.1 ) 


cos ncb+Lj^Ca) sinn4>), 

f(a sin cb^aV cos 6) sin cb = PQ(a) + ^(P (a) 

n^o ^ 

cos n(b+Qj^(a) sin n<b I, 


The coefficients Pj^(a) . . . . , are calculated in the 
usual way. It will be s-ufficient to give the explicit 
expressions : 


( 8 . 2 ) 


1 ^ 

K^Ca) « ^ f (a sin cb , a>^ cos 6 ) cos 6dit) , 

Sn 

= h S f (a sin 6, a^? cos ib ) sin 6d<!). 


Taicing advantage of (8.1 ) we can represent ( 7.7 )j 
( 7 . 8 ) in the following expanded forms; 


^ ^ _ ^K^Ca) - ^^(3C^(a)cos n(Vt46 )+L (a) 
n/o 

sin n(Vt+ <b) ) , 

^ ha F F V" 

if = n(Vt+<t))+Q^(a) 

Sin n (v)t+4)) ). 


Let us Integrate these espresslons in the interval 
t, t ■+ T, within which we consider a, i>, as constant and 
equal to the values a(t), <l)(t). 
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- |rK^(a(t)), 
fePo(a(t)). 

Since T and the increments a(t+T)-a(t), <t> ( t+T)-4> ( t ) 
are small ^ we replace In the left sides hy 

thus arrive at the equations of the first 
approximat ion : 




If we compare with the exact relations ( 8 . 5 ) we 
find that the equations of the first approximation are 
obtained from the exact equations by averaging the righ 1 
hand sides with respect to the time. This process 
duly generalized in the obvious way will be described 
as the averaging principle. 

It need not be said that the preceding reasoning 
cannot pretend to any sort of mathematical rigor. 

For this reason we shall examine in the next chapters 
the questions of the mathematical foundations of the 
averaging principle and likewise the question of 
forming the higher approximations. 

9 . Retumlng to ( 8 . 5 ), if we have a solution in 
a and i and substitute it in (7.1 ), we obtain an 
approximate expression for x. If we choose in place of 
the unknown 4» « Vt+<l> , then ( 8 . 5 ) yields 


dt 


Substituting in ( 8 . 5 ) and (9-1), in place of 


(9.1 ) 
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P thelp expressions from 
gtr 

d a _ ^ f ( (P s in 4 > 


(9.2) 


dt 


2 rr V 


( 8 . 2 ) we obtc-in explicitly 
cos 4> ) cos <i) d<i> ^ 


(9.3 ) 


dvp 

dt 




2'n- 

5 


f (a 


sin 


cos <^ 5 ) sin d 4 > . 


Thus the first approximation to the solution of (6.3) 
will he of the form 


(9.4) X = a sin ip, 

where the amplitude a and the full phase 4 j are to be 
determined from ( 9 . 2 ), (9.5). 

dx 

10. Suppose that F in ( 6.1 ) does not contain ^ 
in which case f will likewise be free from it. Thus 
we will have 

( 1 O.T ) =r f(x). 


and hence instead of ( 9 . 2 ), (9.3): 


( 10 . 2 ) ^ i 


2Tr 


(10-3) ■|^=^(a)= 0 + sin 4> ) sin (i> d<t> 

X 

$ (x) = S f (x)dx 


If we set 


then 

2'n:„, . ,, L l?'^di(a sin ‘i>) ^ 

r f (a sin <i>) cos 4>d<p « — ^ ^ _ o, 

o ^ 

and hence 

da 
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Tlius tlie amplitude of tlie oscillations Is now constant 
a = and so Instead of (10,3 ) we have 

~ do(a)t+e, 

where the phase 6 Is constant and equal to the initial 
value of ^ • 

An approximate solution of (6.3) is then 

(10.4) X = a sin (6o(a)t+6). 

We may say that here the nonlinear character of 
the equation has no other effect In the first approxi- 
mation than to make the frequency depend upon the 
amplitude . 

If we square both sides of (10.3 ) and retain only 
terms in t we obtain 

(10.5) r f(a sin < 1 > ) sin <|>d<i>. 

o 

Since F(x) = ^^^x+tf(x) we have finally: 

P 1 2Tr 

(10.6) CO (a) — F(a sin <j>) 3 ±iii>di>, 

Formula (lo,6) has the considerable advantage that 
the function P enters into it directly and not merely 
through its nonlinear part as it does in (10.5 ). 

11. We will now examine a certain nximber of 
examples . 

( 11 * 1 ) Example 1 . Consider the equation of the 
pendulum reduced for small oscillations (say not exceed 
Ing 30°) to the form 

(11.2) + -£<X- = O. 

dt^ 16 


We have at once from (10.6): 



II. ELEMENTARY THEORY OF THE FIRST APPROXIMATION 


11 


CJ ^(a) = ^ ^ (a aln <t> - tiJlli.) 3in (t>d(t>, 

J. cx. o O 

and 30 approximately 

Cx) ( a ) = 1 - ^ — ) . 


(11.3) 


As the amplitude increases the frequency decreases and 
hence the period increases also. This is likewise 
shown by the approximate formula 


(11 


T . r- 


a 


To take a concrete example for a = 50 ^ we find 

T 


1 . 01 4 X 2fx\[^, 

O 


(11 ‘5) Example 2 . Consider the differential 
equation ( 2 . 2 ) for the oscillations of a shaft. Here: 

ji+jp ’ 

( 11 . 6 ) co^(a) = 


j-j — c(a aln +) aln 4 . (3+ . 

12 O 

To take a concrete case, suppose that M == c(e) is 
represented by the graph of Pig. 6 , or more explicitly 
that 







C(6) 


Pig. 6 
Jh-f-ke, 6>o 
Vti+ke, e<o. 


We find here 
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2 ir 

^ c(a sin 4 ’) 3 in 4 > dll’ = + irka 

o 


and so by (11 . 6 ) : 


( 11 . 7 ) 


2 ,. , '^ 1‘^'^2 . to ^ 


In oraen that this formula be applicable it will be 
clearly necessary that ^ be small. This qiiantity 
measures in a sense the deviation of M from linearity. 

(11 .8) Example 5- Take the case of the electrical,^ 

T2 L -C 

circuit of (2.3) and related equation (2.5)- Assuming 
small, we find by (10.6): 


( n .9) 


2 / _ s A/ . . 5Ba^ N 

o I^A ^ 


from which follows approximately: 

(11.10) ■2^ (a) =\|^1+'2^-). 

12. We will now examine some cases where F con- 
tains ^ 

(12,1 ) Example Consider the equation of 

van der Pol ( 4 , 4 ) where the parameter t is assumed 
small. Comparing with the basic equation (6.3) we 
have here: 


= 


f (X,^) 


^( 1-x 


2 xdx 
^dt-- 


As a consequence we find 

f(a sin aV cos ~ -a(l-a^ sin^ cos 

a^ . o 5 , 

= a(-jjj“ “1 ) cos <!> - ^ cos 3 4 >, 
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and therefore 


2rrv^ 


2it 

i f(a 

o 


sin <i> , av^ cos 4> ) cos <b d4» 


2*^ ' h 


^ f(d sin h, aV cos t ) 3in4> d4 = o. 
o 

Thus referring to ( 9 . 2 ), (9-5)^ (9-^)^ we have in the 
first approximation 


(12.1a) 


X == a sin 4^ 


( 12 . 2 ) 

(12.3) 


da 

dt 




a^ \ 

ir) 


__ 

dt 


From ( 12 . 5 ) we obtain 41 = t +6, where 6 = 41 ^. 
Finally the first approximation is a harmonic 
oscillation 

(I2.il-) X == a sin (t-fS) 


with constant frequency whose amplitude varies in 
accordance with ( 12 . 2 ). By an elementary Integration 
we obtain 


o2^tt 

a^e 


l+ia2( 


e ^'^-1 ) 


and hence finally 

a^e 
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Substituting from (12<.5) in (12.^4-) we obtain the ex- 
plicit approximate expressloxi for x: 

a 

(12.6) X = - 3 ill ( t+6 ) . 

Vl4a|(e«-) 

A trivial solution is x = o wblch corresponds to 
the static regime (without oscillations). It is not 

f 

difficult to show, however, that this regime is not 
stable. Indeed however small the initial amplitude a^ 
may be, it will grow monotonely tending to 2 as a limit. 
Thus the least disturbance will throw the system into 
an oscillation with growing amplitude. 

Prom (12.5) we see also that if — 2 , then 
a = 2 for all t>0. This corresponds to the stationary 
■ regime 

(12.7) X == 2 sin (t+6). 

This "dynamical” regime is strongly stable, for what- 
ever a^(^O), whether large or amall, a(t ) — *2 when 
t — > + 00 , Thus an arbitrary oscillation will tend to 
the stationary oscillation (T2.7). 

The systems of the van der Pol type differ essenti- 
ally from those of the conservative type with equations: 

+ E.f(x) - 0. 

dt^ 

Indeed in the conservative systems as we have seen 
there may occur steady oscillations of arbitrary anipll- 
tude whereas In the van der Pol system steady amplitudes 
are possible only for special values. Physically this 
Is evident from the following considerations : since a 
conservative system neither dissipates nor creates 
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energy, oscillations once started have no reason to die 
down or to grow and so their aniplltudes remain fixed. 

On the contrary in a ” a elf- ex citing” system there is 
creation as well as dissipation of energy and so the 
amplitude may increase if the source of energy provides 
more energy that there is dissipated or conversely. 
There will thus arise a fixed amplitude only if the two 
processes compensate. 

15* ( 1 5 • 1 ) Example 5 . As our next example we 

will taJce Rayleigh’s equation (4.3). Here the function 
f of (6%3 ) will he 

(13.2) f(g) = 

Hence we have 

(15.3) X = a sin (Vt+e), 
as our first approximation, with 

V , e = const . , 

and 

ao -1 

(13.4) ^ ^ f(a^ cos 4>) cos<5 d<t>. 

However we find from ( 15 .^) 

f(a9 coa <|>) «= -aV(A“-|-B^a^>'^ ) cos <|>+^B(av^^) 

C033<|> ^ 

and so from (13.4): 


(13.5) 





20 


INTRODUCT ION TO NON-LINEAR MECHAKICS 


It rollows from ( 13 - 5 ) that the trivial solution 
a — o will be unstable, since A> 0 , and so we have here 
a self-excited, oscillation. Tlie stationary amplitude 
satisfies 


A - 


o. 


which yields 

(13.6) ^ = ^3B- 

Whatever the initial amplitude C=?^o) we have 
from (13.5)' 

a( t ) — f 1 \fhA 

t-^+co T’V 5B ‘ 

Thus whatever the initial conditions the oscilla- 
tion tends to a steady oscillation represented by 

X == (v^t+e). 


ir we desire to learn something not mehely about 
the steady oscillations but about the imtermediary 
regime, we must integrate' (13.5)> which yields 


a(t) = 



Example 6 . As our next example we 
taice an electrical circuit with constant capacity C 
and self-induction L, and containing a non-linear 
element N whose voltage- current characteristic is 




e * P(i), 
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SI 


The differential equation for 1 Is 


(It. 3) 


2 

LC ^ + CP'(I) ^ + 1 = 


df^ dt 

Tliis equation is reduced to the form (6.3) by setting: 




- tf’M ^ ^ ) <^1 


In order to have a clear picture of the degree of 
smallness of the nonlinear element, it is convenient to 
introduce the dimensionless time ' which brings 

the equation to the form 

dT^ VLC dT 


This shows that the application of our results will 
require that the dimensionless quant Ity P ’( i ) remain 

small relatively to tmity* 

We find here that the first approximation ass-umes 
the form 

1 — a sin (v^t+cb), ^ = const.. 


where the amplitude a satisfies the equation 

with 


da 

3t 


aR(a ) 
2L 


2TX 




R(a) =* P*(a aln t) cos^t'dt. 


We see at once from that if R(a) is always 

positive, a(t) — >0 so that the oscillations die down. 

In this case steady oscillations with an amplitude other 
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than zero are ruled out. Referring to (ii^-,5) this will 
certainly occur whenever P*(i))>0 for all i. 

Thus if the characteristic e =* P(l) of the non- 
linear element does not have a falling part (where 
FH1)<0) then the system is dissipative: oscillations 
once started die down. If on the contrary there is a 
falling part in the characteristic then R(a) will be 
positive^ at least for small values of a. In this case 
small amplitudes will increase and small oscillations 
expand, so that the position of equilibrium Is unstable 
and physically impossible, and we are dealing with 
self-excitation. 

Consider the special case where 

e = P(i)=A+Bi + Cl^ + + El'* + Fi^. 

We find then 


R(a) = B + + IbW*. 

We must assume that the coefficient P>o, for otherwise 
beginning with a certain a>a», R(a) will be negative and 
oscillations of amplitude a^>a' will expand to infinity, 
which Is ruled out physically. 

Consider the equation 

B + -^Da^ -f "^Pa^ ** o, 

whose solutions are 


D 

P 


V25P 


zi B 
5 P 


i 
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We will examine the three cases 

(I) B> 0 , D> 0 : (II) B> 0 , IXO: (III) B< 0 . 

In Case I, both roots are imaginar^^, R(a) is always 
positive, hence the system is dissipative. In Case II^, 
if 


fo 


then the system is likewise dissipative. In the 
contrary case there may exist steady oscillations of 
amplitude a^ • However, the system is not self- exciting, 
and oscillations whose initial amplitudes are less than 
a-j dle'‘down. In Case III, the system is self -exciting 
and there is a unique stationary regime for the oscilla- 
tions with amplitude a^a^* 

15. We return to ( 6 . 5 ) and its approximate solu- 
tion X = a sin where a, ip, are given by 


) 

: 1 5 . 2.) 

where 


dt “ 


it “ 


2Tr 


5(a) 

co(a) 


* - 5 f(a sin <t, a.)f cos <t>) cos (t d(b. 


21T 


v>+ Sin d) d*. 


If we square 4.1 and refer to (^.2) we obtain 



2h 
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(15.3) 


2tx: 

P(c. sin 6, 


av' cus 6) sin cb 


and likewise 

if 1 

(15. it-) $(a) = - "grTpS cb)co3 6 d6. 


Thus hy meii-ns of ( 15 - 1 ), ( 15 - 2 ) the functions u), J 
are determined directly in terms of the function P of (6,1). 

16. We will now discuss ( 15 - 1 ) which determines 
the variation of the amplitude in function of the time. 
Ohserve that there must exist no a*^)>0 such that 

l)(a)>0 for a)>a*. 

For if such an a* existed then taking an initial ampll- 

tude , we would obtain in view of (15-T) a(t) — ^+00 

t — *00 which is physically ruled out. 

Referring to (15.1) we see that if the initial 

amplitude a^ is not stationary, i. e., does not satisfy 

5(a) - 0, then with increasing t the amplitude a (t) 

will steadily tend to a stationary determination. 

The tendency of every oscillation to approach a 

steady oscillation points to the special role of steady 

oscillations for all high-frequency oscillatory processes. 

r 

Indeed in such systems the intermediaoTy regime tends 

very rapidly to a stationary regime and hence every 

oscillation may be viewed as practically stationary. 

A noteworthy special case may be mentioned here 

where there are no intermediary regimes, and every 

oscillation is stationary. It will take place for ex- 

dx 

ample whenever the function P does not contain ^ (con- 
servative system). Then (6.1 ) may be written in the 
form 
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( 1 6 . 1 ) 


d^x 

^ + PCX) = 0, 
dt^ 


and 30 direct integration is possible. Indeed if we 
introduce the potential 

X 

U(x) = ^ P(x)dx^ 

o 


then ( 16 . 1 ) yields iraraediately 

( 16 . 2 ) = const. 

Practically^ however, this conservative case never 
occurs and there is always dissipation, hence loss of 
energy, or for that matter there may be self-oscillation 
and production of energy within the system. 

17- We will now consider the stability of the 
stationary oscillations. Let a^ be any root of i>(a) - o 
Then for a very near we will have a = a^ + 6a and so 
from ( 1 5 - T ) : 

= |'(a^ )da. 

This shows that is stable, that is to say, 
corresponds to a stable stationary oscillation if 

(17.1) l'(ap<o, 

while if 


i'(a^ )>o, 

then, the corresponding stationary oscillation will he 
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:..r- 3 t . in particular the static regime (a=0) will be 

' ;jaa *: at 1 e wh er.e r e r 


( - " . 2 ; 0 )]> 0 , 

‘arc 31” this last inequality is the condition for self- 

As we have already seen self- excitation is not nec- 
essary f'tr the existence of self-osc illations in the 
system, that is for the existence of stable stationary 
tscillaticns . For that purpose there must merely exist 
an a . such that ( 1 ) holds . 

An interestlTig case cf frequent occurrence is 
where the system depends upon a parameter An example 
is a series circuit with an Impressed harmonic voltage of 
ar:rlit-:ce u . Inder the circumstances generally 5(a) 
will he a ftnctirn 5:.a, fi). A typical situation is the 
gratn = c of Pig. 7. The dotted arcs represent 

:he -utstahle statior^ry amplitudes, the heavy arcs the 
stalls snes. Through the variation of ^ there may thus 
trise regimes as indicated by the arrows. 



:Fig. T 
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1 c- Having discussed at length and the anpli- 

tvdes we will consider equation ( 15 . 2 ) fcr the frecuencp. 
Corresponding to a given frequency and hence to 0 , given 
period, the frequency is ^ = aj(a), and generally 
f recuency and period depend upon the amplitude. Ihere 
ore, however, important cases in practice when the 
oscillations do not depend upon the amplitude (isochro- 
nous system). An example is when we hao’-e idenoically 

2tr 

\ f(a sin <t), cos <b) sin -t dd = 0 . 


This will occur notably if the initial eqoiaticn is of 
one of two forms 

( 18.1 ) ^ + tf(x}|| = 0 , 

dt 

(18.2) ^ + tP(g) = 0. 

dt*^ 

Notice that ( 18 . 2 ) can he reduced to the form (IS.I) by 

dx 

the change of variable = y. 

Referring to (15) the solution of ( 18 . 1 ) will be in 
the first approximation 

( 18 - 5 ) X = a sin (o^t+d), 

where d is constant and 

(18.4) dt 

2 ^ 

( 18 - 5 ) A(a) == f(a sin t) cos^t dr. 




"3 fi-:'' ecrxatlcr. (5.6) will have a solution 



^9* In the tresent chapter we shall discuss a 
•Tietlti for replacing the first approxiinatlon by one 
«fhich will be stnewhat more accurate, although, of 
crrrse, essentially mere coniplicated* 

Let us examine again the exact equation (8.3) iri 
a, 4 s • ‘The process leading to the first approximation 
ernsisted essentially in replacing the ri^t hand sides 
by their ccr^tant terms. One may thiiik: of these terms 
as cerrespendiTig to slow smooth variation and obtained 
by neglecting the rapidly changing terms represented by 
the trigcnometric functions. In order to account. to 
seme extent for these more rapid changes, we will 
utilize the basic concept of the method of successive 
appTCXiuatioriS . 

Replace then first (8.3) by 



-r^T^TT^-r- 


it 


Ko(a) 

+ ^(K^(a)cos n(v^t+6 

) + 

L,-, (a)slr, rhv 

{ 19 . 

.1 ) 





dt 

dt 

= 

■o(a) . 

+ ^(Pp(a)cos n(\>t+6) 

+ ( 

5 ^(a)sln 

or which 1 

s the 

sarnie by 




da _ 
dt 

oa 

dt 

7^(K^co 3 n(v>t+6) + 

h 

sin n 

( 19 . 

2 ) 






d(t 

dt ” 

dt 

^^(P^cos n(v>t+(i)) + 


sin n ( t +1 ) . 


In these relations a, 6 are the solutions of the 
"smoothed out" equations (8.5)^ that Is to say of the 
first approximations a^ (^, in the sense that we have 
adopted so far. 

Since a, 6 , do not vary very rapidly ^ we will 
integrate the right hand sides of (19.2) as if the a^ 
4 ) In the sioms were constant. We thus obtain 


(19.5) 


— _ n(v^t+(})) - cos n(v^t+(t;) 

a — a "" - ' — - 

nv) 


<t) = 4) 


n(Vt+4)) - cos n(\^t+(|)) 

+ ^ 


The new refined first approximation will be given 

by 


X = a sin (v^t+<t)). 

After some simple calculations, this leads to 







relatlcr-3 Tor a, 4) in tlie 





) uj(a) = \^ + g^Ca) 

ulvaj.em.^y, tc wit:]:".ln qi-iantities of ttie magnitude 

) aj^(a) = + |- (a). 
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20 . Since our results have been obtained "bj me the is 
which have no pretension to rigor, it is necessar^r to 
examine directly the degree to which they satisfy cur 
basic equation (O. 5 ). We find immediately 

2 

(20.1) ^ = -t n(v't+iT)) 

dt 

+ gj^sln n(K*t+(l))i + O(t^), 

where 0(£.^) denotes a quantity of the order of t~ . By 
(19.6): 

(20.2) tf(x,^) = tf(a sin (v^t+i,), 

a.v> cos (v^t+(t>)) + 0(£.^). 


Hence, finally 


(20.3) 


^ = OCtU 

dt^ dt 


±n order words the approximate solution (l9*-4) satisfies 
the initial equation (6,5) to within a quantity of the 
order of 6.^. More explicitly if f(x, x') possesses 
partial derivations of order two or more and if in 
addition the amplitudes are bounded, then 

U-f + + e.f(x,^|<Kt^, o^t<=>, 

dt^ dt 


where K is a constant which depends neither on t nor on t. 
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t::ese : : rjilt 1 :ns ther. ( " 9 . ^ ) satisfies (0,3) to 
1 1 les cf tr.e r-sgnitude of and this 
,nif nnl in fen nii nen"" negative t. 

2 ■ • The euesnicn rf ehe order of magnitude of 
te.e errers rm7 also he treaeec in a different wa^". 

w- first nahe a charge of variables in (6.3) and 
Intmduee new 'nnrnowns a, 5|j , through the relations 


Ein^^ 4. 


} nspeg^sin nt}J 


>-i\ - 

‘V • 


n^-i 


a) oes ^ ■^27 1 "^’i (a) sin 


nV 


.g^ccs n^|J - f^sin n?|J) i . 


*-urer sene sinrle computations it may be shown 
tnat these new variables satisfy the system 

If = w(a) + 

where X, Y are periodic functions of (|j (with period 
Sir}, and regular with respect to e in the nel^horhood 
cf e - o. Notice that (19.6) may he deduced from (21.2) 
by rejectir^g th« terms in 

If we coapare the refined approximation ( 19 . 4 ) with 
the earlier first approximation, we find that the latter 
aerelj represents the first harmonic in the Fourier 
series (19.4). The other hamonlcs will be of the 
order of saagnitude of t. 
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If we examine the stationary oscillations we find 
from (19.^) that they will he periodic with period 

where a is the corresponding stationary amplitude. 
The corresponding frequency is : 

cu(a) « + gip (a) . 

Since the relation between the frequency' and the 
amplitudes is throng the medium of a term proportional 
to we may, if we continue to disregard terms in 
replace in ai(a) the term in t by any other which 
differs from it only by some term in for instance, 
to within terms in we may write 


a>(a) = + ■ g >4 : 


max 


or also 


OL»^(a) 




^max 




22 . We will now apply the preceeding results to 
some examples. 

(22.1) Example 1 . Consider the differential 
equation for a conservative system: 

(22.2) +€.f(x) - 0, 

dt^ 

where f(x) is an odd f -unction. We then verify that In 
(19.5) only the terms remain. As a consequence the 
approximate solution (^ 9 .h) assumes the form 



Tc ::c::-i.i:3A.p rviscKAA^ics 


1. 3 In {oLfZ-rB) 




a) sin n(u^t-4-'6) 


rL>1 


n®-i 


‘^ = JTE 


= W-+|-g^(a), 


;vhene a, 6 are arbitrary constants. 

As a special case suppose that we are dealing witli 

^ V O X 

-f \>^X + tx-" = 0. 
do 

lince d(x) = we rind 

ria sin r) = 5 -a^sln r - -^a^sln 5 ^', 
and hence the refined first approxiamtlon will be 

f 

X = a sin (u^t+e) ~ sin ( 3 i*^t+e), 


where 


iv — + 

C 22 .ti) Sxample 2 > let the basic equation be tbe 
approx inate equation (ll,i ) for the pendulum without 
fraction.. To reduce this relation to the foiTin (22*2) 

we set 


Cas.?) 


f = i'®, = ef(x) 



■X - 
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and ( 22 . 3 ), ( 22 . 4 ) yield here 

X = a sin (ojt+e) + sin 3(iJt+6), 

( 22 . 8 ) 

-\fF<’ - f|)- 

The comparison with the classical series for the 
same quantities shows that we are just ohtainlng t''"e 
first terms of their series, 

( 22 . 9 ) Example 5 . As our next example, we will 
tahie equation (2.2) for the oscillating shaft, he 
choose c( 6 ) as in (11.5), and find as our basic solution 


X = a sin (ccjt+e) + n(a)t +6 ) 

^odd n(n2-i ) 


( 22 . 1 0 ) 


UJ 



+ J 


J1J2 


2 


k(i 




-f 


-rrka^ 


In the present case as it happens, it is not 
difficult to obtain the exact solution, and it is found 
to be 


(22.1 1 ) X - a sin (cut+e) + ^ , 

^ n[(^) V-l] 


(22.1 2) 


= 4 


1 


In this case then the approximate solution may be 
obtained by replacing in the denominators 


n[(^) - ij 



?v1EC-mNTCS 


tr.e frequency a; 'ey its approximation v>. Moreover we 
3 ee thrt (ss.: ; ) yields the accurate expression of the 
first two terms of the exact solution considered as a 

- -a - ^ 

newer series in f~r. 

25. 'A'e will r.DW consider the approximate solution 
{^9.-.} as applied tc a dissipative oscillatory system 


(£5.' ) 

We have here 
(£3.2) 


^ + er(x) ^ = o. 

dt® dt 


»(x,-g) = f(x)^ 


Therefore 

( 25 . 3 ) f(a sin T, cos r) == f (a sin t) aP cos t. 
Introduce the function 

X 

(£3--) P(x) = 5f(x)dx, 

o 

and f om the Fourier series : 

( 25 . 5 ) FCa cos cos ii(t) 

By different iatiag both sides of ( 23 . 5 ) we obtain in 

combtoation with (25.4): 

af(a cos lb) sin i> =*^nP*(a) ain mb. 

If we set *=T + (25.3) yields: 

f(a sin T, ai> cos t) - sin^ ii{^ t+4)+|ff) . 
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Hence (19.-^) yields the approximate solution 

(25.6) X = a sin (»)t+it)) - ^^F*(a)sln n(v/t+ii5+|^) , 

where 4) is an arbitrary phase constant. Here a satisfies 
dt 2 

In particular the stationary amplitudes are the roots cf 

F*(a) *= 0. 

If instead of ({) we introduce another constant phase 
B — i> - then (25.6) takes the form 

(25.7) X « a cos (cat+e) ~ ^^:jP*(a)sin n(Lat+ 6 )- 
where a> ~ 9 . 

Let us apply the argument to van der Pol's equation 
We have then 


f(x) = x"^ - 1, P(x) X, 


and hence 


53, 2 5 

P(a cos i)) ~ - a cos i> — a(^ i )co 3 4 > + ^ 

cos 54 ?, 


so that 


(a) = aC-^ 1 ), F^(a) * ^ , 


P^(a) *= 0 for n 1, 5- 


(23.8) 
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{2^."'} cecomes here 

c 

(53.9) s: = a cca (t+e) - sin 3(t+e) 


where e is ai. arhiirar^” ccnsiant. and a satisfies 


-S' - l§i/ ^ > 

dt ” 2^ ‘ • 

Fcr t?*e statirriar^* oscillations a = 2 ^ and hence 
;23.*c]‘ X — 2 CCS (t-f6) - ^ sin 3 (t+e). 


2^-., He turning to (23.1) we notice that in the 
apprcninaticn cu.der ccnsiceraticn, the frequency is v, 
that is t: say the first term in the expansion of the 
freqtency in powers of t . The second term, the teitn 
in is c in this case* We will now show how to calcu- 
late fcr stationary oscillations the term .in by 
means c f (23."). 

We first observe thiat since a stationary oscilla- 
tion is periodic with a certain period T, we may expand 
the exact solution x in a Fourier series: 

( 24 , 1 ) X = a cos (ti> t-fe)+,c^^A cos n(a»t 4 e)+B sin n(ujt+6), 

n>l ^ ^ 


where a is the amplitude, e the phase of the first 
hamcnlc and ui = On the other hand (25.1) yields : 




dt-b 


+ ef(x)x-^)dt - 
dt 


0 . 


Since x is periodic we ha^e identically 
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T 


S = - Ic^)^dt, ^f(x)xd;-: = c. 

at 


T ^ T 


Hence 

(2^.2) 




Substituting (2^4-. 1) into (24.2) we obtain 


c/(a2 + ^ n^(A? H- B^) = V'^Ca"'^ + ^ (a5 ^ B"), 
n>l ^ ^ n>i *" 




-.2, „a 


and hence 


(24,5) 


(^) 


+^n^(Af + Bb 

n>1 


By compapison oi (25<.7) with (24.1) we see that apppcxi- 
mately 


= 0, B^ = 


-en 


V(n‘=^-i ) 


F„(a). 


We have, therefore, the following expression for the 
frequency of the stationary oscillations: 


(24.4) (^) « 


a n>l n'^-i 


1 + -f-o (-§ — 

n>1 n^-1 ^ 


a • 


It la not difficult to see that (24.4) will hold In 
any case to within tertns of order neglecting 



0 
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tlzei^efore in ( 24 . 4 ) teinis in we obtain tine einrpler 


expression 

( 24 . 5 ) m 


2 


n2-i 


P*.(a) ‘ 
i— ) 


cr finally: 
( 24 . 6 ) Lu 




P^(a) 


2 

) • 


Let 113 apply this formula to van der ?ol*s equation 
( 4 . 4 ). The stationary oscillations correspond to a = 2 
and 30 from (23.8) there follows: 

P*(a) = ^ = = 0 , njf3. 


Hence 

( 24 .?} S5S 1 “ 


IV. cc:L'2TRir<y riQN of the higbbr approximations 

25. In. the preceding chapter we Indicated a way 
to lEsprove the fl 3 ?at approxlmat Ion hy taking into con- 
3 iderat ion the hi gher haimonlca , and we have shown that 
the basic equation (6.5) was satisfied to within terms 
In . Moreover for equations of the form (23.1) we 
gave a method for confuting approximately the stationary 
frequencies to within terms in We shall now consider 

methods for foradng approximate solutions corresponding 
to stationary oscillations whlcdi satisfy (6.3) to within 
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terms in any given power of t. Once and for all we 
will assume that the functions entering in (6,5) have 
all the required derivatives and that the amplitudes 
under consideration are bounded. 

Whenever we state that an expression satisfies the 
differential equation to within terms of we will 
always understand thereby that the error is of order 
uniformly in t for all non-negative t. 

26. Consider first the conservative system 

(26.1 ) ^ + tf(x) 

dt^ 

In this case clearly an "arbitrary” oscillation will be 
stationary. Referring to (22.5) the refined first 
approximation will be 

(26.2) X = a sin (wt+6) + ^ -h^Ca) 

^sln n( 4 ot+ 6 ) + b^cos n(u;t+e) 
n>l - 1 

where a, e are arhltrai>y constants, and gj^ are the 
Fourier coefficients in the expansion 

(26.3) f(a sin t) = ^(hj^sln nr + gj^cos nr). 

Moreover here 


( 26 . 1 ^) 


(a). 
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We will nidify (26.2) 30 as to express it in terms 
oi the ccef ricients of the Fourier expansion 

f{a cos t ) ~ ^ ^ ^ 


F':r this purpose replace in {26.5)TtyT+ thus 

ohtaining 

f ;a cos T ) =H(gn sill n(T + ^) + 11 ^ + §))• 

The ider.tif Icaticn with (26.5) yields 
(2€. = } s^^ = f^sinf?, h^=f^co3f^ 

Substituting these expressions in (26.2) w 6 find: 


X == a sin (ti 3 t+ 6 ) + 




'^o n>l 


fnCOS n (a;t +6 


2 ^ 


r? - 1 


and hence replacing 6 hy the new arbitrary cb ~ 6 - 

we hsFe: 


C26.7) X *= a cos (ajt 44 )) + 


1 - 5 -t 

^2] ^ njl 


f COS n(^ 4 )t+(l)) 


n^ - 1 


This expression is more convenient than (26.2) in that 
It contains only cosine tenns. 

In view of (26.6), formula (26.it-) becomes here 


(26.8) 


(M 


2 \2 



27 , The expressions (26,7), (26.8) suggest the 
following method for ohteining approximations of any 
order. Represent the solution of (26.1) in the form 
X « z(t), where t ^ cot + with 6 an arbitrary con- 
stant and 2(t ) a periodic function of r with period 
2TT. Notice that x = z(r) will satisfy (26.1) if, and 
only if z(t) satisfies the equation 


( 27 - 1 ) 


2 d^z 2 

+ V"^z + tf (z) 0. 


We will endeavor to obtain a solution of (27.1 
such that we have expansions 


( 137 . 2 ) 


z(t) = 2 q(t) +ez^(r) + 


oc 4. 

o ^ 


where the coefficients are to be determined by substi- 
tuting in (27*1 ) and annulling the powers of t. Pur- 
thermo re this is to be done in such a way that the 
are periodic in r with period 2Tr. 

We thus obtain the following recursive relations : 


(27.5) 


Ot JL. 7 

Otjg 4. o2_ 


-f(2o> 


_2 * ^ ^n+r 


P'<Zq, z^) 


-<x ^ 

n+1 


i^ore 


z^y is a polyooailal In 
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Stippcae z,, . . . . , satisfy the 

fli»3-D k + 1 relations of the system. Then clearly 

x = z^+ez,+...+ £^Zj^ 

where 

(a-". 5 ) + ... + 

will satisfy our initial equation (26.1) to withiln 
rems of order t , Tlius x may be considered as 
one required approximation to tbis order. 

The successive determinations of the coefficients 

arbitrary elements which we will utilize 
or remove the secular terms in the solution. Take 


;2T.6) ^ 

as the solution of the first equation (27.5). .The 
second equation yields then 

P ^ ^ 1 

27 . 7 ) ji ( — -i- ) == -f{s. cos x) + oc. a cos x, 

dr^ * ‘ 

Hence In view of (26.5) we have: 




{ 2 "". 8) f coa nT+(cx..,a-f- )coa t 

dT ^ n>2 ^ ' ' 

To avoid secular terms there must be no terms in cos t 
at the ri^it and so we must have 


f,U) 


( 27 . 9 ) 


a 
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Prom this follows for the solution of ( 2 T.S): 


(27.1 0) 



^ n^-i 


In the same way and by an evident induction one 
obtain every 2^ and avoid step by step the presence 
secular terms; the details may be left to the reader 
28. * As an application talre the equation 


(28.1 ) 


dr^ 


4- X 


+ tx 


3 


0 . 


We find here 


(28.2) 


d^z. 


dr^ 


+ 2 . 



(28.5) — a cos r, - T • 

In view of (28.3) the first relation (28.2) beccmes 

= (oc.^a - ^ a^) cos t - ^ cos 3 t. 

*•1 “ 5 t . 

From this follows by the regular application of the 
method : 


dr 2 

Therefore 


( 28 . 4 ) 



JOS + 


cos 5 t, 


thin terms in 


53 (wt+(i5) + (1 - ^ f^) COS 3(u)t+(b) 


+ 1 WoV COS 5(iijt+4)), 


where a, ^ are arbitrary constants and uj Is given by 




The same method may be applied in an obvious way to 


(28, £) ^ + >)^X + tf(x) + (X) + . . . = 0. 

dt^ ^ 


29- Consider again the system ( 26 . 1 ) with f’(x) a 
power series in x: 


(29-1 ) 


f(x) = bgX^ + bjX^ 


Sepe there is no small parameter t . However if -we 
n»rely wish to consider small oscillations then clearly 
f(x) will be small with respect to and inrtherniore it 
consists of a series of terms of increasing oi*der 3 of 
anmll magnitude. This justifies to a certain extent 
the following procedure. Replace ( 26 . 1 ) by 


+ i)^x + pb + P®b,x5 + . . . 

dt® *^2 t 5 


(29-2) 


0 
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where we will consider p as a small paramsuer. Ihis 
equation is now solved as before as a power series ir 
p , after which the parameter p is made equal tc 
thus yielding an approximate solution. 

We will now consider analogous methods for a 
general non- conservative system (6.3). We first moclfv 
our relations (19.^), ( 19 - 5 ), (19-6) for the refined 
first approximation by setting ^;t+6 = and thus 
obtaining : 


(29.5) X = a cos f-F^(a) 


^F^cos r4+G^3ln r4 


n>i 


n - 1 


(29. U ft = I7 Gi(a), f| = co{a), 

(29.5) «(a) = y) + ^ p^(a). 


where F^, are the Fourier coefficients in the 
expansion: 


(29.6) f(a cos Ty -SiP sin t ) « ^(F„(a)cos riT+G^(a)sin nr), 

xi n 


For the stationary oscillations we have 

(29.7) (a) = 0, 

(29.8) 4 » - i 4 ^(a)t+(t). 


where <t> la an arbitrary constant. Thus (29,5) may be 
written more explicitly as 



X 


CCS iaj(c)t 


t < 

V)2 ’ 


-Fo(s 





+ b) 

" 7 ^ 


+ sin ii(^ ( a ) t 
1 


H- 6 ) 
\ . 


Pen censereative s^rstems as we have seen (a) is 
identicallv C and hence the approximate solution (29*9) 
2c.ntains the two arhitrary constants a, <b» We will now 
ctnslder a case where G, (a) is not identically O In any 
Interval zf the variahle a. Assume that G^(a) has only 
simple retts, sc that 11 f'er a certain a: G^(a) = o, 
then the c c rrespond ing G * ( a )^o . Referring to ( 29 * 7 ), 
:'29,9} we see that to every root of G-(a) there corres- 

t ‘ 

pends a certain stationary regime, and that for this 

f 

regime the expression ( 29 -?) depends upon the single 
arhitrary constant 6. 

5C. We will new take up the hi^er approximations 
for non- conservative systems and in the main use the 
same methods as for conservative systems • Write down 
the solution of (6.3) corresponding to stationary oscil- 
lations In the form 

(30.1) X = z (ujt + 4> ) 

where b is an arbitrary constant, uj the frequency, and 
z(x} a periodic function with the period 2-rr. 

We first observe that z(t) must satisfy the 
differential equation 

2 

( 50 . 2 ) to® + ef(z,a>^)=* 0* 

dr 

We newf endeavor to obtain z(-r) and to as power series in 

e: 
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(50.3) 


z(t) = Zo(t) + + . . . 


UJ s= CtJ_ 


■f* ^ + • • m m } 


where is a periodic f-unction with the period 2Tr. 
Proceeding as before by substitution in (50.2) we obtaii 
a recursive system 


'“O 


dz. 




G 


o d- 


° ' dx^ 


where at the ri^t in the nth equation there are only 
terms in , ^n-l their derivatives as well 

as in o^, . . . . , ^n-l * first equation is 

solved as 

( 50 . 5 ) ^ 

where a is as yet an indeterminate constant. Substi- 
tuting in the second relation (50.4) there comes: 

2 d^z^ ^ 

( 50 . 6 ) v> ( — ^ ) = - / (Fy^cos nr + sin nr) 

d*r^ 

+ 2 v>cuj a cos T . 

To avoid secular terms we must have 

F (a) 

G, (a) = 0, CO. = ; — , 

' Svte. 


(30.7) 



T^JTKDDTTnTTnTC TO f^ON-T.l l^An M E O HAI^IGS 


wr-l’^ri ceterzii^ne s- and *o, . Tnis solves (50*6) as 


1 

\>2 


l-F^(a) 


( *c. S } z, = a, cos T 

\ cos nr + G, 

n>i 


sin nr 


n® - 1 


where a, is an indeteminate constant;. Notice in 
particular that in contrast to conservative systems, 
z, is net fully cetemf-ned at the first step. For the 
amplitude a, of the first harmonic will he determined 
ty the condition that he free from secular terms . 

As a consequence of G,*(a)^o and the other assump- 
tion made regardio^ G-j (a) the process may he continued 
indefirh-tely, 

31 . It is to he observed that in the method just 
escpcsed and SLre determined at the same step. 

In ether wczxis, is determined at the same time 

as the function z^Cr). For this reason the joint 
detemir^tlcns should he for 


z^Cuit+h} + • . • + t (^t , 


( 31 . 1 }■ 


V + 


^ P N+1 

+ C. LU- 


N+1 


that la to say with x up to the order and a> to the 


order t 


N+1 


For instance if N — O then 


3: *a a cos' (ctft + <b), 

G,(a) = o 

^K&iXctx Is our first aRproxlmation, obtained, previously h^ 
tihe ai^raging ixpinoipl©. For N -« 1 we have 
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( 31 . 2 ) 


+ 


X = (a+ta^ ) cos (<.»t+ 4 )) + -^f-PQCa) 

^P^(a) cos n(-.t+<b) + Gj^(a) sin n(a>t+ 4 ;) 
n>i - 1 


^ = i> + I^F^(a) + 


By and large, the result 3 are the same as bercre 
as rar as the first step of approximation goes, except 
that they are now obtained systematically and net by 
some special device. 

52. We will now consider the same questions not 
merely for statlonaiy oscillations but also for the 
general case, that is to say, for oscillations which 
need not be stationary. 

The form of the refined first approximation sug- 
gests looking for a solution of (6.5) of the form 

(32.1) X z(<ji, a), 

where z( 4 »,a) is a periodic function of (|i with period 
2 Tt, and where 

(32.2) « A(a), — ^(a). 

By differentiating and substituting in (6.3) we obtain 

-2 ^ a2 + Is a 

3*1 3+3a 3^t aa 

( 32 . 3 ) 

If 35^ +. <>^2 + "+ If A) * o 

It la clear that If we find, z. A, «», satisfying (32.3) 
to within any particular oirfder of magnitude In t then 
(32.1), provided that ( 32 . 2 ) holds, will satisfy (6.3) 
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tc witliln zne same erder* To find the required expres- 
sions cf 0, A, ui, we will set 


z(4}, a) = z^Hi, a) a) + . . . , 

{^s.k)<^ A(a) = tA^(a) + t%(a) + . . . . 

I ui(a) = + €. ( a. ) 4- • * • • • 


where we as some periodic in + with the period 

gtr. These expressions are sucstituted in (32.3) and 
yield the system 


(32.5) 




+ z. 




+ 


Z, 


-f(Z, 






1 a^paa 


2 vnj 




The first is solved as 

( 52 . 6 ) 2^ = a cos 4 >. 

fe could equally start with any other solution, for 
instance = a sin but this would hot introduce any 
essential change anywhere. Substituting then in the 
second equation of (32.5) we obtain: 

+ z^ = 2 i^{n^a cos + + sin 

-I<-n COS aln nif>). 


(3S.7) ( 
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To avoid secular terms, we must have 

(32.8) =P^(a), 2(>A^ = G^(a) ^ A, 

and the resulting equation for Is solved as : 


(32.9) Zi = -^l-FoCa) 


^ F^(a) cos n 4 f + 

y-s I 


n>l 


" 1 


The process continues in the obvious way. We thus 
obtain in succession 








up to any index n. For Instance if we have reached the 
value n ~ N, then we obtain a solution to within the 

TJo. T 

order t of the form 

(52.10) X = a cos 4) + €.z^(tp,a) + . . . + ,a) , 


where a, i|i , satisfy 


(32.11 ) 


^ + e.^Ag(a) + 


+ t^Ajj(a), 

. . + t^/\(a) . 


Thus for N * l we obtain precisely the formulas for the 
refined first approximation. 

35. Referring pow to the first equation (32.11) we 
see that the stationary amplitudes are given to within 
order by the solutions of the equation: 

(53.1) 0. 
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Let be a root of G^(a). Under our aasimiptlons 
it is not a double root, and so (53.1 ) may be solved 
in the form 


a +a^t + * • • ♦ 

In the first approximation we will have 

G>(a^) 

The stationary regime under consideration will be stable 
if 


~jG| (a) + 6.^A^(a) + • . . + t^A^(a)<Co, 

and unstable otherwise. Since clearly the left hand 
side is 


+ o(e.^) 

the question of stability for t small enou^ will depend 
u|»n the sign, of that is to say we obtain the 

same criterion as for the first approximation. 

Cfeaerve also that the eqxiationa for the Nth 
approximation like those of the first, show that the 
aiplltx^e a will Increase or decrease monotonely ap- 
proaching fpota a3x>ve or from below the nearest stationary 
a^littiie according to the sign of ^ for t =* o. 

In geneml, one must eit^jhaalze the fact that, 
axcept for certain singular cases, the relations for the 
ftMt approxlBatlon provide the same qualitative Indi- 
catlona for the starting of self- oscillations as the 
hJgJber ^pprooclaations. Generally speaking, the higher 
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approximations provide quantitative rather than new 
qualitative information. In view of this and cf the 
difficulty of confuting the higher approximtian, it 
usually quite sufficient to obtain the first apprcxir. 
at ion. 


V. LIHEARIZATIQH 


In the present chapter, we will first of al 
endeavor to obtain suitable interpretations for the 
equations of the first approximation. We begin by 
writing the basic differential equation in the form 

m ^ + kx + tf (x,^) = 0, 
dt*^ dt 

where m, k, are positive. This system has two well 
known interpretations, the one mechanical, the other 
electrical. 

We have obtained as the first approximation a 
solution 

(5i,2) X = a cos 


where a, ^ satisfy 

r 

I tia t 

2Tri>'in 


2tl 

\ f(a cos (b. 


sin <t)) sin 4> di 



« a»(a) 


( 3 ^- 5 ) 



Il^TRODUCTI-N 


rc meckaitics 


21^ 

(5*.. 4 ) = ‘'^ + -ili I c|>,-a^3in t) 

COS <t> d6. 

It Is to te kept in mind also tliat the first approx Im- 
■^tion (54.2) represents the fundamerital. harmoni c in the 
expression of the refined first approximation (see for 
Instance formula (29.3)) which satisfy ( 3 ^. 1 ) to within 
Gi^er t , 

Let us introduce the functions of the amplitude 
k(a), A (a) defined by 

2'*^ 

(5t.5) A = ^ f(a cos 6, sin 4 )) sin <t) d6 


, 2^ 

(3a. 6) k = k + sincb) cos 6 d(j>. 


In terms of these quantities the equations ( 5 ^- 5 ) for the 
first approximation take the form 


( 3 ^- 7 ) 


da 

dt 


A 

m 


a. 


Lt 



As a consequence, we obtain by a direct If lengthy 

coinputation: 

m + A — + Ex = 0 (e.®) • 

dt^ dt 

He isay then say that the first approximation (5^-2) 

2 

under consideration satisfies to within the order t 
the linjpar equation 
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2 

( 5 '+. 8) m + A ^ + kx 

dt 

In short In the first approx im5.t ion the o sc illations 
of the non-linear system londer consideration c.re ecoi-/- 
alent to those of a linear system with a dissipation 
coefficient A and a spring constant k. The approxina- 
tion is to the order t , that is to say neglecting 
qnantlties of the same order as when we formed the 
first approximation. For this reason we will call a 
the equivalent dissipation coefficient, and k the ecniv 
alent spring constant* The linear system (3t.S) will 
also he said to be equiva-lent to the assigned system* 
From the comparison of (34.8) with the given equation 
(34.1 ) we see that the former arises from the ncn-lineo 
system by replacing the non-linear term or restoring 
force of the mechanical analogy 

(3i^.9) F = 

by the linear term 

(3i^.10) k,x + A 

where k^j — k - k. 

Let ns remark also that ^ = 2m disslpatic^^ 

decrement in the equivalent linear circuit, and 
the proper period of its oscillations, to* within the 
order of f ^ . 

33. We may conclude then that the equations (34*7 
of the first approximation may be derived as follows: 
Linearize the system by substituting for the restoring 
P of ( 34 . 9 ) the restoring force P^ of (54.10) where a, 
k^ are defined by: 





2'rT 

J 5= - — ^ \ : 'a :;is i, -T.y sin 4^) sin 6 d6 

•' T! T k' ■;:: 

2 rr 

f 

' :^r.S; \ I'a 2 :'S r), -aP sin 6) cos (b do. 

The eq*nrl"'n 2 for 5^ as demanded by a linear 
syster 


-< 53 , 


^ ^ <u; 

dr 


and obey ^.re precisely rrnse of the first approximation. 

Tbs f:r:rrl process ^rst described will be referred 
n: :.s o^lirr rrlnolple of linearisation. 

J-. Tnao is the physical significance of linear- 
Izorlcn’’ To answer the qrestion we will have to have re- 
cr'Tse :: -,n electrical system. 

he first recall certain concepts familiar In 
electrical engineer Ir^. Let 

^ e(t) — S cos ost, l(t) = I cos (t^t - oc) 

be a hamcnlc o’cltage and harmonic current In a given 
dm: It. The angle ^ Is the phase- lag of l(t), and 
cos is the power-factor of the system. The complex 
represent at Ives of e, 1 are 

;i. 2 ) nt) = Ee'5“'^, ■T(t) = j =\pr. 

Dent tir.g temporarily by x the conjugate of any quantity 
Xf if T » is the period of the oscillations, then 

T 

iS^(t) ir (t)dt = p.-jp^ 

^ cA X 


13 C. 3 ) 
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where , the power of the 
effective power) is hr. owe es th^ 
es the reective ocwer. '".e e-iso 


a vs t e.v , ' ve ' s ^ 




■ 3-.0 


O = hS e(t)l(t)it 

= if i e(-)i(t - ^).ai 


These lest expressions rmp senoe to def 1: 
o.ny periodic e, i. 

Now (34.1 ) represents the ir.cticn cf 
subjected to the force -kx-e.f. Assvning 
h--moriic end of period T, the ne'an power 
ncti’/e power will he 



( - 5 ) 


T 

m ^ ( hX-?-t I ( X , X ^ ))X ^ Ot 

c 


By an obvious o-nalogy we may introduce here also a 
reactive power 

T 

( 36 . 6 ) = 1 5 (toc+£P(x, x'))x'(t - pdr. 


If we impose upon the linear system ( 3 4. 6 ) the conditi 
that its active and reactive powers be P_j to with 

p ^ 

terms in t , we obtain precisely the values given by 
(35.1 ),(35.2). 

37 . Another physical Interpretation may also be 
obtained quite directly as follows. Substitute the 
harmonic oscillation x = a cos (i^t+e) In the relat loros 
(34,9), (54-. 10 ). For this harmonic oscillation the 
equivalent linear force F.j will be likewise hani»nlc 
with frequency v) - Let 6.,, denote the phase and 
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IT o- 


= I- CCS (v>t-fci)^ ) . 

ar f'rce will 'oe periodic but with various 
cse free;: ere ies will be multiples ol v^. 
?.:.err = l r.arrr.eric be I cos (v^t+b). If we 
rplltuce and phase of and of this funda- 


I, = I, 6, - 6 

then we ictiiin relations which yield again (55.1), 

In point cf fact in expanded form the linear 
:rr?e will le 



>,a CCS (i^t+e) - v^Aa sin (v^t+6) 
wl.ile one f '.ndanental harmonic of the non-linear force 


— f . a CCS T, -av^ sin r) cos t dr j cos (v^t+6) 


2 

r 

- 7 f ^ f (a CCS T, -aM sin r) sinrdTl sin (v^t+ 6 ). 

If we eouECe the two it is but a step to ( 55 . i ) , (55 . 2 ) . 
The prccess described for obtaining k, A will be 

referred tc as the principle of harmonic balance. There 
is nc difficulty in showing that it is in fact equivalent 
tc the first procedure for deriving (55.1 ), (55-2 ) . 

^t la irnportant to observe that there is no reason 
whatever to derive the differential equation for the 
oscillations before linearizing, the system. Indeed, in 
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many cases (especially for more or less compli? „T:ec 
oscillatory systems) it may actually ce rrare ccnver.ien' 
to dispense with the formation of the differencial 
equation, or to form it only afterwards and cc linei-rire 
the system directly from the data. The casic face is 
that we are dealing with systems which do net differ 
too much from harmonic systems. 

38^ We will now consider a few examples. 

(58.1) Example 1 . Suppose that we have a particle 
subjected to a non-linear spring whose effect is de- 
scribed by F — f(x). Then for a harmonic cscillaticn 

X ~ a cos the fundamental harmonic in F will he 

2lT 

^ f(a cos t) cost drl cos (v^t-re). 

Therefore by the principle of harmonic balance v;e may 
replace the non-linear spring by a linear spring whtse 
spring constant is 

2fT 

k(a) = ^ I f(a cos ct^ ) cos 6 do . 
o 

(38.2) Example 2 . Consider a circuit with an 
iron core and let i be the flux and current w'ith 

(38.3) I = i’(i) 

as the relation between them. If the current is .tar- 
monic : 

(38.4) i = a cos (vJt+e), 


then the fundamental harmonic of the flux will oe 
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2tT 



CCS d6 I cos (i^t+6). 


Ttievefore^ 'hj the principle ol" harmonic balance we' may 
replace (I 0 . 5 } hy the equivalent linear relation 
$ « L^i where 

2it 


By analogy with linear circuits, we will call the 
equivalent coefficient of self-induction. 

(55-5) Exouople 3 - Suppose that we have an elec- 
trical series circuit with the same inductor as in the 
preceding example and In addition a linear inductor 
with self-induction coefficient L and a capacity C. 

By linearization we obtain an equivalent system with 
coefficient of self-induction L + and capacity C. 
Therefore the frequency is approximately 

1 1 , , V 

ay (1 ) w 

V(L+L^)C \/lC 2L 

C3S,6) Example h . Consider an electrical circuit 
with a non-linear element N and characteristic relation 
e = -P(i), If the current .is again given by (38.t) 
then the fundamental harmonic of the voltage will be 

2TT 

{ 38 . 7 ) I - I P(a cos 6) cos (J) d6 ] cos (v>t+e), 

o 

and so the non-linear element N may be replaced by a 
linear ^Ims^nt with characteristic relation 
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srt 

(3E.8) e = -Rgij S * --• 


This assuiTies, of course, tbo-t tne circuit is suclt tr.ut 

the oscillutlors ure nearly harmonic. 

If is nosltive then the circuit acts as an 
e 

ohmic resistance and abscrhs mean power to the .amt ant 


of 



If on the contrary is negative then the circuit 
generates power to the same amount in ah sc lute value. 
The system is then said to have the characteristic cf 
a generator. 


VI. APPLICATION OF S7?^0LIC MEOKCZS T: 


39 . Let us introduce the linear operator J wtatse 
domain are the sines and cosines, and which is defined 
"by 


j sin 4 i)t = cos wt, j cos wt — - sin 
so that in particular 

(39.1 ) f — 1 • 

It l 3 a consequence of (39-1) that • j has the character- 
istic values ±1 (1 = V-l ) and related characteristic 
functions e“^“*. Moreover If iCz) is a ftmctlon 

of the complex variable z and J(i) — A+IB, A and B real, 
then $( j) _= A+jB. Furthermore If A+IB = re^“ then 
J(j) = re^t We also have 







COS + i + 


T 3 Z-¥t 


:3 . s cos ‘^-ot 

rio-3- If fCt) Is hanoonic and of 


:: r.en 


^f(t) 


:,eri‘ 



1 

to ^ • 


• Traslder now a linear conductor to whose 
lerr'lails la :.pplied a hamcnic (sinusoidal) voltage of 
fre::-e:u2i Kirchrff *3 law will yield a linear differ- 

erui^l er,.:.tirn with ccnst: 3 nt coefficients for the 
?-rrem I't,. In view of ( 59 . 1 ) if there is a harmonic 
2: then it will satisfy a relation 


s the impedance of the 
ve for the complex current 

Z ( J m)T — “’e . 

It is ::ften ccnverJLent to introduce also the Inverse 

cperattr 

(- 2 . 5 ) A(ju>) = 2^3137 



kncwn as the admittance of the conductor, hfore generally 



let 2Z "be a network with two terniir^ls and let e, 1 
the sane meaning. Klrchofl-s laws yield then the 
similar relations and sc 21 has an Impedance ant an 
admittance. 

As a simple example if an Indnctcr L,a tap at 
C, and a resistance R are ccrnected in series^the 
impedance is 

(1.0.3) Z = Lji*. ^ H ^ ^ 

while if they are in parallel the admittance Is 


(10. i.) 



The concepts of admittance and impedance which 
proved so Important In the theory of alternating cn 
rents have been extended in recent years to ether 
branches of physics notably tc mechanical and actus 
systems . 

Consider for instance the motion of a part id 
governed by the equation 

(1I-0.5) m + icx = fCt), 

or with the velocity v = as unknown, by an equati 

-i- 

(ko.6) m^+Av+kS vdt = f(t). 


Consider on the other hand an electrical circuit go 
emed by the relation 

t 

L H + Ri + H “ 


( 40 . 7 ) 


e 
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where L, R, C, e, have 
view cf vhe ccnplete f? 
differential eqvatlcns 
Pierre Curie (see his 
the fcllcwing analogies 

V - -.V- ^ C 2 0 il 1 at icns 

Risplacenent x 
Velocity v 
Fence f 
Maas rr. 

Friction Coefficient A 
Cpring Constant l-c 


(4o.7) one may, with 
p. 16ky 1891) establish 

Electrical Oac illations 

Electrical Charge q 
Current i 
Voltage e 
Self-induction L 
Resistance R 

Inverse of the Capacity C 


the usual interpretation. In 
•rmal identity between the two 

( 40 . 6 ) 

brks” 


Chis is in the main the electrical-mechanical analogy 
utilized in nedern acoustics. 

In connect ion with mechanical systems, one has 
frequent occasion to consider rotating systems. The 
oasic differ-ential equation for the possible oscillations 
of such a system with one degree of freedom will be of 
the fem 

,*^0.3) J f 4- n — rce = M. 

dt^ dt 

It will be seen that it is obtained from ( 14 - 0 . 5 ) if 
X is replaced by the angular variable 6, the velocity by 
the angular velocity, the force f by the torque M, the 
mss by the moment of inertia J, the friction coeffi- 
cient A by the friction moment referred to the -unit of 
velocity n, and finally the spring constant by the 
coefficient of hardness c. We thus have the following 
analogy between retatlng mechanical systems and 
electrical syst^is; 
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Rotating: Oscillations 

Angular Displacement 6 
dQ 

Angular Velocity ^ 

Torque M 

Moment of Inertia J 
Braice- torque n 
Hardness c 

Consider now the harmonic oscillations of these various 
systems . 

Pop the electrical system (io. 7 ) we have 

(^0.9) = Zi, Z - Liuj + R + ^ . 

On the other hand hy the operational method applied tc 
(^ 0 . 6 ) we obtain 

(40.10) f = zv, z = mjuj + A + ^ • 

By analogy we introduce the mechanical impedance z and 
in association with it the mechanical admittance T == y* 
41 . Consider a linear electrical network without 
impressed voltages ^ If we examine the possible exist- 
ence of harmonic self -oscillations of a given frequency 
(Aj, Kirchoff *3 laws yield a linear homogeneous algebraic 
system of equations with a determinant A(juj) rational 
in Joj, A necessary condition will then be 

(‘H .1 ) a( ji") = 0. 

The roots of ( ii-i . i ) will specify the acceptable f re- 
quencies . 

\4i. 2) Example 1 . The net X consists of a single 


Electrical OsclllatlGns 

Electrical Charge q 
Current 1 
Voltage e 
Seif-inducticn h 
Reslst'Unce R 
Inverse Capacity 7 ^ 



- 
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closed circuit with characters L, R, C. Then (41,1) 

cecccies 


(41 . 3 ) 


Z( ju>) == Lja> + R + 


O, 


whichx is equivalent to 


(41 . 4 ) 


R 


0 , ^ 


Vlc 


Thus harmonic self -oscillation is only possible when 
the resistance R = o, and then uj has the value indi- 
cated. These are, of course, well known facts. 

(4-!, 5) Example 2 . Let the net 21 have two term- 
tmls and let there be impressed a harmonic voltage 

e — E cos (u)t-fot) 


at the terminals . The corresponding complex voltage is 

■g = 


and ac the ccinplex ciirrent i is defined by 

Z( = e . 

If we have Z(iuj} = !Z(3a>)i e^f^ then we find 

(41.6) 1? = , 

!Z( ja»)i 

and hence 


(^1.7) 


i - I cos (wt+«*-p), I = I . 
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One recognises liere the well known relations cf alter- 
nating current theory, 

(^‘S.S) Example 3 . consists of two circuits with 
characters (L, R, C), (L^ , ) and coefficient of 

mutual induction M. If T, are the (complex ) c'urrents 
then Kirchof f ^ s laws yield 

= 0 , 

(i+1.9) 1 _ 

(L^ jo) + - M>i = C, 


and so 

(in. 10) A(juj) = ■*■ “ “• 


Setting 


(in .1 1 ) 




q = 




P p 

the roots u)^ of (4l.lo) (considered as a quadratic 
in ) are given by 


(41 .12) 


gUl)2i\|(v/-2-wf)2^4q2g2\)2 

2(1 - q^) 


Pi-ovided that q^ i . The admissahle frequencies will 
be to^ , tjL)^ if they are real • 

42. By combining the preceding developments with 
linearization, their range of application may profitably 
be extended to nets with non-linear elements, as we 
shall now show. 

Returning to the first exan 5 >le the system may under- 
go an oscillation of the form 
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( 42 




1 = A.e 


-6t 


CCS («U3t4-6)^ 


where 
( *+2 . 2 ) 


6 = 


R 


2L 


uO) 


j^\l’ - ^ ■ 

\[lC V 1+L 


R 

A 3 S‘r:rlr^ row H small and considering ^ as of the 
first order of smallness we will have as a first approz- 
imaticn = "7=^ soid thus nearly harmonic oscillations 
of f requency^^St* Tc neutralize dissipation let us in- 
sert in the circuit a ncn-linear element N with char- 
acteristic e = -FCi), such that the "instantaneous” 
resistance Fh^i) is of the order of R and not always 
positive in the range under consideration. Referring 
to liS-S) we replace N hy an equivalent linear element 
wltit character!"^ tic e ~ -R^i, where 

2tt 

(42.5) ^ 6)cos 4) d6 

2Tr 

= ~ ^ F*Ca cos 6)sin^ 6 d4> 

c 


where a is the amplitude of i. Clearly R^ is of the 
same order as P*(i) and hence of the same order as R. 
Ihe equivalent linear system is a circuit with char- 
acteristics (L, G) and so this time as in (41.2): 

(42.4) -4^ =s O, R + R 

oiC ^ e 


Re(a) 


Vlc 


(^2.5) 


~R, CO 



VI, APPLICATION OF SYMBOLIC MSTHODS TO LllTEPMZZPfriZ:: 


Its occurrence is clecrly impossible imiless is 

not always positive, i, e., 'unless P(i) Has "rolling” 
parts . 

Consider more generally a linear net 51 with in- 
pedaxLce Z( jco) short-circuited on the same non- linear 
element N as above. This time linearization yields 

( 1 2 . 6 ) Z ( j co) + = 0 . 

Hence if Z( jto) = X(aj) + jY(ai) we will have for a sta- 
tionary oscillation i = a ccs (ojt+e): 

(42.7) Y(c^) - 0 

(t2,8) "" -X(co). 

The first relation determines co, and then the second 
the omplitude a. 

If the characteristic of the non-linear element N 
were of the form i = f(e) we would proceed similarly 
with impedance replaced by admittance, and resistance 
by conductance. 

The operator method duly generalized may be 
applied to non- stationary oscillations. C-enerally 
speaking in a linear system a non- stationary oscillatic 
is of the exponential-harmonic type: 

(45*1) X = Ae”*^^ cos (u>t4-<|t). 

It satisfies the relation 

(i^5.2) ^5^ = (-<5+»f(t). 


Hence all oui* arguments may be extended to exponential- 
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harmonic oscillations prc^=^ideG that is replaced 
everywhere hy the operator p = -6 + j^- For example if 
a network has for characteristic equation (4i .1 ) then 
its non- stationary (understood exponential-harmonic) 
oscillaticna are governed cy 

(t3-5) a(p) = 0 . 

Thus for the same network JSI as before with non-linear 
elen^nt N we will have 


ik3»k) Z(p) + R^(a) = 0, p = 

Having determined by this relation 6 and o) as functions 
of a, the elements of (^ 3 .'S) will be given by the equa- 
tions of the first approximation 

(^3.5) H = -da, uj, 4 , = a,t+(b. 

For in the first approximation the solution assumes the 
form (13-1 ) with a = = urt+cj) and this implies 

{^3.5). 

Suppose in particular that SI consists merely of the 
series circuit (L, R, C). Then Z(p) = Lp+R-t^ and 
(43.4) reads 


Hence here 

( 43 . 6 ) 


Lp + (R+Rg) + ^ 0- 


R+Rg{a) 

<’<^5 2L ^ 


There will be self- excitation in the system if (4-3.1 ) 
does not die dcwn when a Is neai zero, i. e. If 4(oKo 
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T1 


or if R^(o)<R , R as -R. Thus R* (here ~R) Is a critical 
equivalent resistance for N such that helow it the sjstets 
is self-oscillatory, above it is not. 


VII. MULTIPLY PERIODIC SYSTEMS 

kh. Up to the present the oscillations under con- 
sideration have been taken so to speak one at a tlsie. 

If the system is linear and several frequencies are 
admissible, say if the characteristic equation (ll.l) 
has the roots , . . . , then there are possible 

stationary oscillations 

{kk,^) -^h “ 

and so by the principle of superposition (for linear 
systems ) there is a stationary solution 

{kk,2) i — a-jCos (u)^t+<b) + . . . + a^cos • 

It is fairly clear that this principle may not be applied 
to a linear system equivalent to a given non-linear 
system. Under certain conditions (reasonable smallness 
of suitable parameters) some progress may still be made, 
as we shall now show. For simplicity we limit the dis- 
cussion to the case of two oscillations. Two distinct 
situations will arise according to the presence or 
absence of resonance. 

^5. Consider then a linear system 51 with two ter- 
minals across which there is connected a non-linear 
element N whose characteristic we write as 
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e == ~F(1) = -ef(i) 


wl'iere t will serve to gauge the deviation from linearity, 
'e will suppose also that f is a polynomial. In prac~ 
ti':re this is a rather mild restriction if f is contin- 
uous 1 the usual case), since it may then be arbitrarily 
'und uniformlxr approximated by a polynomial over any 
closed interval* 


iuet 


z:>) 


oe impedance of the network and 


1 0 ^ 


imental (natural) frequencies. We will then 


Lve 


Where g! Jau ^ ^ o, h = l, 2. 

Ass'uming new t small let the non-linear system 
aam.io the cscillaoions represented in the first approxi- 


Nt-i) = a^cos (a;t-f(i)^), 

where are arbitrary constants and where o 3 ? 

n 

- to the order , Set now 

■U5.u) 1 = i.+ig = a^cos (a>t+43^ ) + a^cos (ojt+d)^). 

The resulting voltage in N is 

(*^5-5) e = -C.f{i) =£ -€.f(a^cos ) + a^cos ) ) • 

Consider now the double Fourier expansion with respect 

to 4i^ , d> 2 : 




(45.6) = -tA^^CGS ( (ina;, }t-r.T4 , -r.t ' 

applied to tlie linear network of in:pedance Tnis 

is where we must distinguish between two pcssiile slou- 
aticns . 

16. NON H^SCNAIvZ? 5Y3ZZ2^ . Suppose flrso that nans 
of the expressions 

(16.1 ) 22 + (n+'! )ct;_ 

except those corresponding to m = 1 j, n — 0, and 
n = 1 are of order at least t. This will cental 
if none of the frequencies 

( 16 . 2 ) (m+i)ctj^^ + nto^Q, -f {n±: 

other than those corresponding to (m, n) ~ 
are zero. Referring then to (li*5) it is seen that 
under the circumstances for (m,n) ^ (1,0), (0^1), 
induces In 51 a current 

^A^cos ( imco^ 4-n«g 

|2;(3(maij+n£^g) f 

Z( +na>g) ) ** [Zle'^mn, 




( 16 . 5 ) 
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^ --jr* '± VJ 


Th, e .!iirip 11 1 ud e 


i 


mn 


is 


i Z ( j ( mc*3^ -hULd^ ) 1 


Under the nss'Limptlon +nuj2dirf ers by a finite quantity 
f ren ^mn order t. 

On zhe ccntrnry for instance Z( jco^ ) la of order t 
snd 3 0 similarly Iq^ y is finite. Tlius to tbe 

order t the current generated will Toe 

(:,d.5) 1 = 1,0+ Iqi = (a.^ti+6T-fiTo) 

+ loiCOa ). 


Uhls muso he the same as (45.4) to within the order t. 
Ihe idenoif ication yields first ^ = Pq,, = O to the 
crder , a condition already satisfied. Then we must 


I,, Q, 


-2 == ^01 


-e.A 


1 o 


-tA 


01 


lZ(jcu^)\ 


! Z(jui^) j 


?rcm this follows readily that if we set 

2-it2Tt 




'eh 


1 


2t,2a^ 


^ ^ P(a^co3 <b^+agCos (t)2)cos (bj^dcb, 
( 11 = 1 , 2 ) 


o o 


then 


( 46 . 8 ) 


® -^^el 
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Thus in the first approximation the non-linear o' 
istic e = “^f(i) = -tfCi^+i^) may ce replaced ly 
linear characteristic (46 •8), We may therefoc'e 
pret the replacement of the non-linear element IT 
equivalent element with the characteristic y.:sr ■ 
as a linearization of the system when there is n 
resonance . 

As in the case of a single oscillation we 
write down here also a refined first apprcximti 
i, representing it to within the order and i 

he, neglecting no i^ term: 


(46*9) i = ^^xnn "" (oj^t+ci)^ )+a2C09 (ca^ 


|Zj(m<4j+na)2)l 


(46.10) 


mn 




47 . RESONAM^ SYSTEM . Suppose now that one 
frequencies (46.1 ) is zero, or uj^q ~ ^^ 10 ^ where 
are relatively prime and -- is one of the fracoio 
It is clear that we may assume wioho 
striction co^u)^ , or r^s. This time we will have 
finite Fourier s-um: 

f(a,^G03 (cL)^t+6) + agCoa (w^t+is^)) 

=.y]i^co3 

If we set 

®in = 


( 1^7 - 2 ) 
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then we have 


The current 


(47.1.) 


e = -^X%co3 (f^tH- Sj^). 

Induced hy e^^ in is 
i = -^^nC03 

1 r^/ V , 


(^T.5) 


Z( 


The amplitude of i^^is 


(47.6) 


^ = 


iZ( jfco^ )! 


and it is finite for m == r, a, of order 6 otherwise. 
Thus 


i = ir+ig = IgOoa (cu^ t+e^- )+i^co3 (oj^t+e^-p^,) . 

The identification with i = i-j +^3 yields here 

^ 3 ^ ^2 ~ ^r' ^3~^s ~ ^ ®2?'”^r ~ ^ 2 * 

The linearization aaauraes this time the form 

e -(Z^i^+Zgig) 

mfmrG are (coni^lex) impedances whose computation 

offers no jmrticnlar difficulty. 

The ’•refined” first approximation for i^ 

iMition to the order ie 


or appi»ox- 



VIII. 


rNFLUETvCS OF PERICIII - 

(47.7) i* = cos (o3^t+6^') -f a^cos t+d^) 

(f^,t^e-p^) _ 

n¥r,s IZ(j|w^)| 

In addition to the harmonics uj^ , it will ccnraln 
others of the form > which may be of smaller fre- 
quency than uj^ . Thus we have here so-called derr-:lti- 
plication of frequency, a property of considerable 
p rac t i cal import anc e . 


VIII. INFLUSNCE CF PZRICIIC lI5-Ir.3.e:rZS ■ 


48. Up to the present we have concentrated upon 
isolated systems, not subjected to any exterior dis- 
turbances . As an example of a non-isolated system we 
will discuss the equation 

2 

(48.1) m ^ + kx = 

dt 

where k, m are positive, t is small and 

(48.2) f(t,x,-g) = fo(x,^) +^(±-*(x,||)co3 A^t 

4. A^t), 

where the sum is finite and f^, f*, fj* are jxslynoiiiials . 

The mechanical interpretation of (48.1) is obvious. 
As equivalent electrical system we may choose a series 
circuit with current 1 - Inductor m, cai»clty g. 



condenser ci3a.rge x, and non-linear element N whose 
characteristic is 


Cii-8.5) 


e - tr(t,x,i). 


Since for t small the system is quasi -harmonic we 
shall apply the general concept of linearization. For 
£ « 0 we may choose 


(^ 8 . 4 ) 


(48.5) 


X = a sin^ (t^Qt+(i>), 

1 = ^ = auj^coa (uj^t+i). 


( 48 . 6 ) 


o Mm 


For t small hut we will consider the above formulas 
as approximations and substitute them in (48.5). We 
have then 

(48.7) e *ef(t, a sin (uj t+<t>), auj cos (uj t+6)). 


Since f^, f*, t** are polynomiala we have finite 

Fourier a rms : 

(48.8) f^Ca sin ^^aiu^cos ) «^Jf^(a)co 3 k:»|*+gj^(a) 3 in jk 

aiid alallarly for f*, f*. with f* , f*g, and g*g 

as the coef f loi«ita . Hence 

(48.9) e -XcfjcCoa k(u)^jt+<l))+g^alii k(uj t+(|)))+ I. . . . 


the wwrltteaa terms conta ining sines arid cosines of the 

ranges of k and of the A 





VIII. IliFLIlEtirCE OF PERIODIC DISTI 3 RBAMf!Ef?. 

finite. Here again we nmat distinguish between reson- 
ance and non- resonance accordingly as k:ui^+A_ does 

or does not hold for some (k,A^). 

49 . NON-KESONAM* SYSl^EM . 11113 is the case where 
no frequency is itself, or where (k-i 

whatever k, in their ranges. Then the only haim)nic 
of frequency uu^ in e is 

(1^9.1) e, = £(f^(a)coa (ui^t+cfc) + (a)3iii (uj^t+i)). 

In view of (48.6): = Z^J?, where 

(49.2) Zg = ;^(f^ (a)-Jg^ (a) ). 

o 

By the basic principle governing linearization we re- 
place the non-linear element N by an equivalent linear 
element with characteristic 


The equivalent linear system has then the Impedance 

characteristic equation Is z(p) 

» Z^, or explicitly (see 45 ): 

© 

“P + i ” (a))^P = -<5+j“- 

Therefore : 

in{-6+ujj) + = il^(f,(a)-3gi(a)), 

and this yields to within the order t 

^ ** 2m^a 


( 49 . 3 ) 



32 


If we coirblne wit^- the equations of the first approxima- 
ticn C see ^ 3 ) * 


dt 


-i a. 


dt 


ss UJ 


jr 


we rray replace the latter by 


da __ t ^ 
dt ~ 2 mu^ ^ 1 


(a). 




— Lu - — - — (a) = ui(a). 

dt o 2ma^a ^ \ 

4> = uj^t+4). 


The related first approximation for x is . 
(49.3) X = a sin it . 

Introduce now the expressions 

(^ 9 . 6 ) 


A = (a), k| = -| Si (a). 


auJL 


Since we have to within the order 

o 


uP 


and 30 Ci9-^) inay be replaced by 


(% 9 . 7 ) 



Wixm (^9*7) we deduce to within the 
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order the relation 

(^9-S) — “ + ^ 4. (jc + k^)x « c* 

dt" dt ® 

Notice that (^4-9, 8 ) depends solely upon the tern cf f\ 
Since 


the linearization and associated first approxlirjatlcn 

may he obtained by applying the averaging nrccess to f 

dx 

(averaging as to t, as if x, ^ were independent vari- 
ables), and replacing f by the resulting f'anctlon f * 

To sum up then: as regards the first apprexin-- 
bion and related linearization we may replace f by 
fo(x,^). Since this last fiinction does net contain t 
explicitly, we have a situation already considered* We 
merely recall these properties : 

(^9 > 9 ) The stationary amplitudes =51^0 are the solu- 
tions of A (a) = 0 . IT a„ is such a solution then the 

o 

corresponding oscillation is stable whenever 

and unstable otherwise* This assumes, of coursfe, that 

the derivative at a^* 

(if9.lo) Self- excitation occurs when and only when 

A (o)<o. 

(50) By way of exaiople let us apply the preceding 
results to van der Pol’ a equation with a forced oscilla- 
tion: 

( 50.1 ) 

where as uaxial £)>o. To reduce (50.1 ) to the form (% 8 * 1 ) 
set 
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8A 


( 50 . 2 ) 


+ b sin OLt, b = 


E 


1 -ot^ 


Then x satisfies 

» ^ o fix 

( 50 ^ 3 ) + X « eCi-Cx + b sin a't)'=)( — + bo. cos ott) 


Here 


m 




- X 


2 \ dx 

^ at 


^2 - k 2 

A = (1-7;; — ) - 


Hence the first approxiinatlon is 
X *« a sin (t-»-<t)) 

( 50 .l^) <t> = const. , % (1 - f- 


Therefore there is self- excitation when and only when 
b ^<2 and there Is a stable stationary amplitude a = 

V 4-2b^. The corresponding stationary solution of 

( 50 . 1 ) is 

( 50 . 5 ) y aa b sin ott + \j 4-2b^ sin (t+<t>). 

¥qt b^ 2 , X =* 0 is stable and, so 

( 50 * 6 ) y « b sin oit 

la a stable forced oscillation for ( 50^1 ). 




51 - RESONAOT SYSTEM. To siaiplifT matters we will 

dx 

suppose f of the form f , where f(T,u,^7) is 
periodic in r and of period 2fC. The basic equal Icr. la 
then 

2 

( 51 . 1 ) m ^ + kx = £f(o(.t,x,^). 

dt*" dt 


We suppose now that 

LU =r ^ oo + tCl 
O 3 

where is an irreducible fraction. As usual we set 

( 51 . 2 ) X = a sin (~od:+<t>) 

and replace P = t f by the equivalent linear force 
= “kgX-A^ . By identifying the fundamental hamcrJ.cs of 

tf(cxt,a sin (|at+4)) 7 ^ cos (foft+6)) 

and 

-k^a sin (^at+c|)) cos (§at4-4)) 

we obtain to within the order 

_ 2tr „ 

^ r f( 3 T~ — 4i, a sin rr^auucos rr) sin rrdr 

■ e tra J ^ r^ ^ 

2tr 

* ( fCsT- a sin i*T,aui,cos rr) cos rrdr, 

*ftauj_ 0 ' s ' o 


The equivalent linear syat^ thus 
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(51.5) m + A— + (k + k )x = o. 

dt^ dt ® 


Thus here to within the order f 


k + k^ ]£ 

(51.!.) a = ^, ^ u.^:i + ^), 

and the equations or the first approximation are (see 

13 ): 


(51 .5) 


^ - -iA p 
dt 2m^ 




dt 




Ir p 

since £ « ujt . 
m o 

(51.6) It is to be observed that if f(T,u,v) is a 
finite trigonometric sum of terms sin kr, cos kr, then 
•unless (r,s) is in a certain very limited range the t 
term in f (t,x^-^) does not influence k^, A and as re- 
gards the first appi^oximation it may be suppressed. In 
that case we ai*e back to a previous case where f Is of 
the form fCx,"^:^-) (no periodic disturbance). Hourly 
speaking it means that the resonances that count occur 
within a limited range of values ( r , 3 ) . 

Motice in particular that if ^ then the fre- 

quaacy will be very near u)^ and hence very near Thus 
the application of a dlstiirbance of frequency oc may 
induce 'an effect of frequency This Is known as sub- 
harminic resonance or deinultipllcatlon, and the property 
has been extensively applied especially in radio tech- 
nlcpie. 

Jjn point of fact, iK>t only will the frequency — 
appror but also certain isiltlples which depend upon the 
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nature of tine polynomial f. 

(51.7) Let U3 apply tlie preceding ccnslderatlnn 
to (50.1 ), the van der Pol equation with hamcnlc dis- 
turbance. For s>3 we find the same situation as In !5C 
and nothing is changed. For s « 5 we must replace the 
second relation of ( 50 . 3 ) by 



Here again self-excitation in x arises only for 
with a limiting stationary amplitude 

a = -2b + \jh + 2b . 

The corresponding stationary solution of (^c.i) Is 
given in the first (not refined) approximation by 

y = (-2b-N4+2b^)sin (^+4>) + b sin od:. 

Notice that when b^ == ~1— (weak disturbance) then 

i-or ^ 

the subharmonic dominates the harmonic . 

For b)>2 the situation is as In ( 50 ) and there is n 
subharmonic . 


IX. CQMPLE 14 ENrS 

52. Wb will first discuss a somewhat different 
manner of obtaining the hi^er approximations from the 
p2?ocedure indicated in Chapter IV. Consider then the 
differential equation 

(52.1) d!| ^ X, t) 

^ fit 


mdiere for t sufficiently simll we have a pcwer aeries 
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I'epres ent at ion 

(5S.2) f(t, X, t) =21&^j^(t, X, 

in idiich ia a polynomial in x, sin t, coa t. 

Regarding uJ we aasume explicitly that it la a positive 

irrational number. 

Introduce now new variables a, 6 defined by the 

relations 

( 52 . 5 ) X « a sin e, ^ = auj cos 6. 

This enables us to replace (52.1) by the system 


{32. h) 


da 

dt 



€.f(t, a sin 6^ au cos B, t) cos 6 

- f(t, a sin B, ato cos B, t) sin e. 


Ilnder our aasuniptiona we also have 

/f^(t, a sin ao) cos e) coa e = F(a) 


(52 


• 5f( 


—X 

m +n. ^0 




f^(t, a sin B, a cos e) sin e « 5(a) 


V 


m^+n®^0 


*C„e 

nm 


3 (m 6 +nt ) 


ifeere the smm are finite > Iiet now 



IX 


CCMPLEMETITS 


ii(a 


(52.6) 


Y“ -•(r:e+r.t] 
i“'~’ ~ X—'^nn. i(vr.uH-r.) 

^ \ -■(ir.e+nt) 

v(a, e, t) NU •. - , , 


We verify at once the relations 


(52.7) 


= f o ® ■ P(a) 

It ^‘^11 = 


We inti*oduce now in place of a, 6 new variatles a. ^ 6 , 
defined by 

a.^ u(a.j , 6 ., , t ) 

“ ujIt ^^^1 " - li Wja, , e. , t: 

By substituting the expressions ( 52 . 8 ) for a, 6 in 
( 52 .^ 4 -) we find 

/da^ t ^u(a^, e^, t) da^ >a(a, , t) de, 
at“ + w 35“ ' Se, 3T^ 

"^uCa, , e. , t) c , „ 

^ “ & f(t,a 3 in e, aui cos e, )co3e 

(52.9)< 

3e, t "SwCa^ ,e, ,t) da, "^wCa^^e^jt) de, 

' ’ 5 t- - 65 ar- + Te^ dt 

^ _ lo- f(t,a sin e,au> cos e, t) 

aln e. 



e = 



m p, w 


5C 




:z mechanics 


llrr*ln£uting 6 tj meaxis of (52*8) axud solving for 

G6 ^ 

we octain after some simplifications: 


1 "" t ^ . t, t), 

(52.1 cK 

^ = “J- $(a^) + t^'sCa^, t, t) , 


wnere fcr t very small we have expansions : 



S - 




t ), 


wlGh 5-v polynomiaa.s in cos 6^ , sin 6^ , cos t, sin t. 


R-, = ?. (a. ) 


3„ = i,(a, ) 

where the s'mis are finite. The same reasoning may now 
he repeated with f^ cos B, f^ sin B replaced by S^, 

etc. The final result may he described as follows. 

For each n there may be written a system of differential 

equatio.n3 in a^, e^: 


rr,2 

m +n 


m^+n^^o 


ej(me^+nt ) 

inn 

I 

M^e j > , 


dt 

(52-n ), 

d©, 

lit 


-F(a^) + . 4 €X-i(^) 

+ . . . . 

A / o \ ^ 4- c \ 
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In these equations have the same prcperti- 

as R, S. Pupthermore = R, = 3 and IT 

r(^) =Xe®R^^^, 

then 

C"' ’ - "kOh.) 

4''"’ -'Sk<»k> 

We may also assimie that If t is so small that terms of 
order of may he neglected then we have the follow- 

ing relations for the nth approximation: 

X = V 

where 

u; L- §(a^) + + 

The method just described for obtai ning the suc- 
cessive approximations is very direct and lends Itself 
rather well to an estimation of the error consequent 
upon neglecting certain terms ^ 

Another obsei*vation to be made is that the systeta 
(52.12)^ may be deduced from (52.11 ^ be replacing 

tbe latter by the constant term in its expression as a 

double Fourier aeries, th^ rejecting tenra of order 
^n+l 




n 


(52.12) 


dt 


^ de. 
dt 


n 



::rrHcc:rri:!; to ncx-liksar rwTSCHANics 

A Iasi reriark regarding th.e process just described, 
is tiat ccritinued Indefinitely It does yield formal scl- 
as power series in t , but -unfort innately as shown 
oincare, ihe series are generally divergent. Thus 
■ cannot be utilized directly to investigate the 
structural properties of the solutions. 

Passing new to an entirely different type of 
censideratiens we will discuss the following problem: ~ 
'Ahat indieaticris do the approximations provide regarding 
the ?xaat solutions? 

Tabling first J52.12), dropping the index l : we 
ha'je the system 


( ’ ) 


= F(a) 


St = ■ JJ5 


The corresponding first approximation (the earlier 
^refined” first approximation) is 

X = a sin 6 -f~ju(a,e,t)sin e - v(a,e,t )co 3 ei . 


If we set 

(35*3) r^(t,a sin cos e) ^ 

then we find readily in place of ( 55 . 2 ) 

(a) 

X = a sin e +& X 


(53. 


uj^ - 

(n® + (m^-1 


Suppose now that a*^io Is a sin^ile solution, of 
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(55-^) P(a) = 0. 

Then and we will suppose sxpllolol;-' 

(??.6) P'Ca*)<o. 


Tne other case vrauld he dealt: with hy replacla,: every- 
where i: by -t. 

It follows from (53-5) that (1 5.0 rias the s: lo- 
tion 

(53-7) a = a*, 6 = t/t+4,, t)=uj - ^ Via.*;, 


where i> is an arbitrary constant. Ihe ctrresttnding 
stationary solution given by (55-3) Is entlicltly: 


(55.8) X = a*sin (v^t+jl^) +t. 


Z -P f 




u> - ( muH-n ; 


Thus it is of the form 

(55-9) X = z(t, v^t) 

where z(6, is a continuous periodic function of 
4) with period 2 Tr in each and depends upon b. In parti- 
cular z(t,Pt) will be quasi-perlodic for all irratior-al 

v). 

In view of (*53-6) we .see that every solution of 
(12) for which the Initial value a is near enoit^ to 
a* will tend with t — * +«to one of the stationary solu- 
tions of (55-T)- 

Thus we may assert that approxlimte solution 

(53.4) whose Initial values x, ^ are neai* enou^ to the 
initial deteimilixations. of the aK> 3 ?oxiJiiate stationary 




solution will ter^i with t- + oo to one of these station- 

f 

ary regimes. 

One may prove the following result : The property 
just formulated for approximate solutions (representa- 
tion by quasi-periodic functions of the form (53 -9) 
and properties of stability) belong also to the exact 
■^elution of the differential equation (52.1 at least 
whenever t is sufficiently small. 

This important property shows that the investiga- 
tion of any particular approximation (for Insfance the 
first) obtained by the methods which we have repeatedly 
discussed, has a meaning not merely for purposes of 
approximation but nay serve likewise to give heuristic 
indications regarding the structural qualities of the 
exact solutions. 

The proof of this theorem has been given at length 
in Memoire No. 16 of the Bibliography. We will merely 
discuss here two special cases which will serve as a 
strong indication regarding the nature of the theorem. 

Consider first the case where w is not an 
:^nteger and (53.5) has the soliition 0 with 

F»(0)<0. 


We have thus = 0 for m =3^ o. Hence the corres- 

ponding stationary solution (55-8) assumes the form 


(54.2) 


X 


« E 



u? - ji= ' 


where 

< 54 . 3 > 


u I t At, 0, at. 

O 





It is inunediately evident that this solution, which Is 
independent of the constant of integration, is periodic 
with period 2tr. From the physical point of view It 
corresponds to forced vibrations. 

In view of (5^.1) the approximate solution ( 5 ^. 2 ) 
is stable, and to be precise: an arbitrary approxima- 
tion ( 55 • ^ ) whose initial values x, are sufficiently 
smalls, will tend to the approximate stationary solution 
i^k.2) for small enough t and with Indefinitely increas- 
ing t . 

We will now establish the same property for the 
exact solution of (52.1). For this pui-pose we obser^/e 

first of all that for given iniulal values x , x« the 

o c 

solution of (5^*2) may be represented as a power series 
in t . We will then have ; 


X(t) = (XqC 03 Uit + -^ Sin un;) + tX(t,x 




jx'it) = -(Xj^coa u* - x^ui sin wt) + t) , 

where X(t, x^, x^^, t) ia an analytical fxraction regular 
for aufflclently artiall t. It la clear that will 

be periodic with period 2n, if, and only if, we have: 

(5^.5) x( 2 'n:) ~ = 0 , x*( 2 tr) - x^ = O. 

From ( 5 ^- 5 ) we obtain the following relations for x^. 




\x^(co3 2iru»-i ) + -^sin 2frtti + e X{2Tr,XQ,x^, t) — o 


^x^LU sin 2traH-x^{cos 2fri«>-t > -i- eX£C 2 ti:^x^,X^, t) *• O. 
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For ^ = 0 these r*el-.ticns r*ave the trivial solutlorx 
= 0, xi = with a ner.~ zero jacchlan 

i(ct 3 2Tiar:), ^ air 5tTU3 

p p 

( 54 .’*’) : = (cos 2 rtuJ ~1 ) +sin 2 iTa>^ 0 . 

; -14' sin 2 iiaj, ( cos 2 Cuj -1 ) j 

Frcn this we !oay conclude that ( 5 ^. 6 ) has an analytical 
solnticn for t S'nff iciently snail. Snhstltuting this 
solnticn in we ohtain an anal^vtical expression 

for the periodic sclntion of the differential equation 
( 52 .^). Evidently the constant term in the expansion of 
this periodic solution is equal to 0 . An elementary 
computation yields for the next term the expression 

If we continue with the same reasoning which is 
used in the well known method of Poincare - Liapounoff, 
we will readily see "that the periodic solution under 
consideration is stable. For the characteristic expon- 
ents are in fact FUc)±j owing to and since 

t. is always assmed positive, their real part is nega- 
tive. This proves stability. 

55 . The second case which we shall now treat is 
where the equation (55.5) has a non- zero root which 
satisfies (53-6} and where furthermore f does not con- 
tain explicitly the variable t, that Is to say where 

(55.1) fCt, t) If) . 

In this case the appi*oximate solution (53.^) assumes 

the form: 




( 1 Ho^ 


gjine 


( 55 - 2 ) 


X 


a sin 0 + 
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where 


2rr 

(55-3) fnj(a-) = f^Ca sin 6, a.w cos e)e~-"'^ae. 

Thus the approximate stationary solution is of the for: 


(55.4) X = a*sin (\^t+6) +t 


f (a* ) * / - 1 *. . ^ 




( 1 )iiJ^ 


where 6 is an arbitral^ constant of integratior 


(55-5) iJ = LU - $(a*). 

It is a ready consequence of (55-0 that for indef- 
initely increasing t, every approximate solution (5:. 5) 
with initial values near enough to those cf (55.“^) 
tends to one of the approximate solutions. 

It is also immediately clear from (55*^) that in 
the case under consideration the approximate staticnai»:%^ 

-I 2'JX 

solution will be periodic with a certain period ana 

physically speaking corresponds to free ncn-dissipatir-g 
oscillations . 

To establish analagous properties for the exact 
solutions, we observe first that since here the functions 
R, S, do not contain explicitly the variable t, (52.10) 
may be written in the form 

da _ £F(a) E.^a)R(a.e, £) ^ 

J • de ^2 _ |^(a) + t^u>S(a,e,t) 


Since 


P(a*> == o, f'(a«Ko 
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it f::llcws fpcm the theorem of Poincare -* Liapounoff 
that (55.^) has a periodic solution 

C55«T) a ==TT(6 j t }j TTca> o) = a*, 

where rT(e, t.) is an analytical function of t , regular in 
the vicinity of o, and with period 2Trwith respect to 6. 
On the other hand we have from (52.10) 

h§ ^ ^ iia) , 4. t^S(a, 6, t) 


and hence: 

!w“''+eTT(e, t)! de = at, FT* 


- i 

uaCuj- + t^S) 


By integrating we find 

(55.8) (w - t: + 5 , S - const. 

where 

2n 

M^(t) =^ 117(6 

If we set 

« 55 - 9 > Z ^Te^M^(e) 

then hy the ingjlicit function theorem, the solution of 
(55*8) Biay he put In the foim: 


(55.10) 


o - {^t + 4.) + ♦, t) 
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wiiei>e e ) Is an ana.lytical function^, regular Ter 

sirfriciently small and with period 8tt with respect t 

6 . 

Thus in view of* foxroulas (52.3), ( 52 . 8 ), (55.^), 
(55.'^0), we may conclude that in the case under censid'” 
eration there is an analytical periodic solution: 


X = z(Mt + ij, t), 

where la an arbitrary constant, U and z(6, £) analyti 
cal functions regular for t sufficiently small, and in 
addition z is periodic in 6 with period jf z and 

are expanded in power series and terms of order 
hl^er than one neglected, we obtain again ( 55 .I), 

( 55 . 5 ). 

If we write down the variation equation corres- 
ponding to the periodic solution ( 55 . 10 ), it is easily 
seen that one of the characteristic exponents is c, 
while the first term in the expansion of the real part 
of the other la P ^ ( a* ) , hence for t small enou^A the 
real part is negative. By reference to the theories of 
Poincare - IdLapoianof f , we see then that any solution of 
the differential equation ( 52 . 1 ) whose initial values 
are sufficiently near to those of the solution ( 55 .’ 0 ) 
may be represented in the form 

(55*11) X « z (i> t + , ce t ) 

where z(e, h, t) is an analytical function of h near 
h o, whea?e furtheimore z(e, o, t) equal z(e, 8 ) and 
finally p la a characrerlstio exponent. As for f , c, 
they ane constants of Integration, with c sufficiently 
small. 

It Is thus clear that for indefinitely Increasing 
t, the general solution (55*ii) tends to the periodic 
solution { 55 * 10 ). 
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EEMTA 


p. 1, 1. 11, Jg instead of 0^, 0^ 

9 , 1 . 15 should read : 

enas of the form tx (a trigonoaaetric function). In the 
p. 25, equation {16.2) should read: 

(16.2). +U(x)« const. 

p. ^5, equation (28.2), delete the subscript o from the 

p. 51, equation {32.3 )> at end of equation insert = 0. 

p. 33, eq’uaticn (32.8), G^(a) instead of 

p. 56, equation (3^*7), k instead of k 

p. 59 « equation {36.14-), at end of equation insert dt. 

p 60, 1. 19> k instead of k 

p. 65, equation CI4O.5) should read: 

2 

(2*0.5) m ^ + A + kx - f(t), 

dt^ 

p. 75? equation (%5-6) should read: 

{%5*6) f(a^ cos cos ig) = J cosCmj)^ + latg). 

p* 77 ? equation (146.9)? the first m should have sub- 
script 1 

p* 85? 1. 13 should read: 

a aiii(~ott + 4 >)? cosCf-oct + <l>)) 
a as 
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